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PREFACE 


This  investigation  was  conducted  by  the  US  Army  Engineer  Water¬ 
ways  Experiment  Station  (WES),  Hydraulics  Laboratory  (HL),  under  the 
direction  of  Messrs.  F.  A.  Herrmann,  Jr.,  Chief  of  the  HL;  H.  B. 

Simmons,  former  Chief  of  the  HL;  and  J.  L.  Grace,  Jr.,  Chief  of  the 
Hydraulic  Structures  Division.  The  effort  was  supported  by  the  En¬ 
vironmental  and  Water  Quality  Operational  Studies  (EWQOS)  Program,  Work 
Unit  VIIIA.3  (CWIS  No.  3 1 604 ) ,  entitled  "Evaluate  Alternatives  for 
Aeration/Oxygenation  of  Hydropower  Releases."  The  EWQOS  Program  is 
sponsored  by  the  Office,  Chief  of  Engineers  (OCE),  US  Army,  and  is 
assigned  to  the  Environmental  Laboratory  (EL),  WES.  The  OCE  Technical 
Monitors  were  Mr.  Earl  Eiker,  Dr.  John  Bushman,  and  Mr.  James  L. 
Gottesman.  Dr.  J.  L.  Mahloch,  EL,  was  the  WES  Program  Manager  of  EWQOS. 

The  study  was  conducted  under  the  direct  supervision  of 
Mr.  Jeffery  P.  Holland,  Chief  of  the  Reservoir  Water  Quality  Branch 
(RWQB);  Dr.  Dennis  R.  Smith,  former  Chief  of  the  RWQB;  and  with  the 
coordination  of  the  US  Army  Engineer  District  (USAED),  Savannah,  and  the 
Tennessee  Valley  Authority.  Messrs.  Steven  C.  Wilhelms,  Michael  L. 
Schneider,  and  Stacy  E.  Howington  prepared  this  report.  Assisting  in 
the  testing  were  the  authors  as  well  as  Messrs.  Holland,  Charles  H. 

Tate,  Jr.,  Hubert  R.  Smith,  and  the  Clarks  Hill  Reservoir  personnel. 
Messrs.  Gary  Mauldin  and  James  Gallagher,  USAED,  Savannah,  provided 
guidance,  technical  assistance,  and  support  during  the  field  studies. 

The  report  was  edited  by  Ms.  Jessica  S.  Ruff  of  the  WES  Information 
Products  Division. 

COL  Allen  F.  Grum,  USA,  was  the  previous  Director  of  WES. 

COL  Dwayne  G.  Lee,  CE,  is  the  present  Commander  and  Director. 

Dr.  Robert  W.  Whalin  is  Technical  Director. 

This  report  should  be  cited  as  follows: 

Wilhelms,  S.  C.,  Schneider,  M.  L.,  and  Howington,  S.  E.  1987. 
"Improvement  of  Hydropower  Release  Dissolved  Oxygen  with  Turbine 
Venting,"  Technical  Report  E-87-3,  US  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  Miss. 
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IMPROVEMENT  OF  HYDROPOWER  RELEASE  DISSOLVED 
OXYGEN  WITH  TURBINE  VENTING 


PART  I:  INTRODUCTION 


Background 

1 .  Hydroelectric  power  generation  has  proven  to  be  one  of  the 
most  attractive  energy  sources  available.  Past  energy  shortages  have 
spurred  growth  in  the  number  of  hydroelectric  facilities,  and  numerous 
existing  and  proposed  sites  for  hydropower  projects  are  being  evaluated 
and  developed.  Hydropower  is  presently  meeting  about  12  percent  of  our 
nation's  energy  needs.  Additionally,  several  thousand  potential  hydro- 
power  sites  have  been  identified,  and  the  attractive  attributes  of 
hydropower  have  resulted  in  evaluation,  design,  or  construction  of 
hydropower  facilities  for  many  of  these  sites. 

2.  From  the  standpoint  of  energy  resource  conservation,  hydro- 
power  is  very  attractive  because  the  energy  source  (water  held  in  the 
reservoir)  is  renewable.  Hydropower  is  also  flexible  from  an  opera¬ 
tional  standpoint.  Changes  in  power  demand  due  to  daily  peaking  and 
seasonal  fluctuation  dictate  the  need  for  a  rapidly  responding  energy 
source.  Hydropower  generation  can  usually  be  stopped,  started,  or 
changed  in  a  matter  of  minutes  by  simply  controlling  the  flow  rate  of 
water  through  the  turbine.  This  provides  nearly  optimum  compatibility 
with  peaking  demand.  If  the  supply  of  water  is  abundant,  hydropower  can 
also  be  operated  continuously  to  meet  baseload  power  demand. 

3.  Hydropower  is  considered  one  of  the  cleanest  major  sources  of 
electrical  energy.  However,  adverse  environmental  impacts  resulting 
from  a  proposed  or  modified  hydropower  project  must  be  evaluated  and 
techniques  that  minimize  or  mitigate  damage  to  the  environment  must  be 
developed . 
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Problems  and  Concerns 


4.  A  frequently  cited  problem  associated  with  a  proposed  or 
existing  hydropower  project  is  the  release  of  water  with  a  low  dissolved 
oxygen  (DO)  concentration.  This  problem  is  typically  the  result  of  low- 
level  releases  from  a  density-stratified  pool  coupled  with  an  in-lake 
oxygen  demand.  Due  to  the  heating  of  surface  waters,  a  reservoir  strat¬ 
ifies  such  that  a  surface  layer  of  warm  water,  the  epilimnion,  resides 
above  a  layer  of  cooler  water,  the  hypolimnion.  DO  concentrations  in 
the  epilimnion  are  generally  high  due  to  the  extensive  transfer  of 
oxygen  at  the  air/water  interface.  The  hypolimnion,  may  however,  become 
oxygen  deficient.  The  presence  of  density  stratification  acts  to  in¬ 
hibit  vertical  mixing,  thereby  limiting  the  transfer  of  oxygen  into  the 
hypolimnion.  When  a  hypolimnetic  oxygen  demand  is  coupled  with  this 
absence  of  oxygen  replenishment,  deterioration  of  hypolimnetic  water 
quality,  often  to  the  point  of  anoxia,  occurs. 

5.  Several  problems  may  develop  under  anoxic  conditions,  such  as 
the  dissolution  of  trace  metals,  release  of  nutrients,  formation  of 
hydrogen  sulfide,  and  depression  of  pH.  Hydropower  intakes  are  often 
located  in  the  hypolimnion,  resulting  in  poor  water  quality  releases 
downstream  during  power  generation.  Depending  upon  the  severity  of  the 
DO  deficiency  in  the  release,  it  may  be  necessary  to  employ  one  or  more 
techniques  to  enhance  DO  concentration  in  hydropower  releases. 

6.  The  retrofit  of  an  existing  flood  control  or  other  nonpower 
project  with  hydropower  has  produced  a  number  of  water  quality  concerns 
(Wilhelms  1983).  At  many  nonpower  projects,  significant  reaeration 
(Wilhelms  and  Smith  1981)  (often  to  near  saturation)  occurs  in  the  high- 
velocity  regions  of  open-channel  flow  through  the  outlet  works  and 
stilling  basin.  The  incorporation  of  a  downstream  turbine  and  pres¬ 
surized  conduit  results  in  the  loss  of  this  reaeration.  While  the 
impacts  of  this  loss  of  reaeration  are  often  site-specific,  a  change  in 
release  quality  from  highly  oxygenated  to  near-anoxia  (as  described  in 
paragraphs  4  and  5)  would  severely  impact  the  downstream  environment. 

7.  A  number  of  Corp  of  Engineers  (CE)  projects  have  been  designed 
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with  multilevel  intakes  in  order  to  withdraw  water  from  different  levels 
in  a  stratified  reservoir.  This  provides  a  means  of  releasing  water 
with  various  temperatures,  DO  concentrations,  and  levels  of  suspended 
sediments.  The  addition  of  downstream  flow  control  resulting  from 
hydropower  retrofitting  may  limit  or  negate  this  capability.  Many 
potential  concerns  regarding  water  quality  have  been  encountered  and 
identified  at  existing,  proposed,  or  add-on  hydropower  projects.  Based 
on  these  problems,  guidance  is  needed  for  the  design  and  operation  of 
hydropower  reaeration  techniques  to  enhance  release  water  quality. 

Study  Objectives  and  Scope 

8.  A  wide  variety  of  techniques  are  available  to  improve  the  DO 
of  water  released  from  hydropower  projects.  The  system  most  effective 
at  a  specific  hydropower  site  will  depend  upon  many  factors,  including 
the  degree  of  DO  enhancement  required,  rate  of  release,  turbine  type  and 
operation,  upstream  and  downstream  water  quality  objectives,  and  avail¬ 
ability  of  economic  resources.  Selection  of  the  "best"  system  must 
involve  weighing  the  costs  and  benefits  of  each  technique  with  regard  to 
site-specific  concerns.  The  initial  objective  of  this  study  was  to 
investigate  means  of  enhancing  hydropower  releases.  This  investigation 
led  to  extended  evaluation  of  one  technique,  turbine  venting,  in  terms 
of  dissolved  gas  uptake  and  the  costs  incurred  by  altering  turbine 
operating  characteristics.  The  results  of  this  study  should  provide 
general  guidance  for  the  applicability  of  turbine  venting  to  projects 
similar  to  the  study  site  at  Clarks  Hill  Reservior,  Georgia.  However, 
caution  should  be  exercised  in  application  to  other  type  turbines,  low- 
head,  or  small  projects.  As  a  backdrop  to  the  documentation  on  turbine 
venting,  an  overview  of  techniques  for  improving  reservoir  releases  is 
presented . 


PART  II:  POTENTIAL  TECHNIQUES  FOR  RELEASE  ENHANCEMENT 


9.  As  mentioned,  the  most  frequently  cited  adverse  impact  for 
proposed  or  existing  hydropower  projects  is  the  release  of  water  with  a 
relatively  low  DO  content.  A  wide  range  of  techniques  are  available  to 
address  the  problem  of  low  DO  concentrations  in  project  releases.  These 
potential  solutions  vary  greatly  in  terms  of  economic  impact,  opera¬ 
tional  complexity,  and  degree  of  influence.  The  following  discussion 
outlines  the  many  techniques  which  have  been  identified  to  improve  DO 
releases  from  reservoirs  and  lakes.  Additional  details  on  many  of  these 
techniques  were  presented  by  Bohac  et  al.  (1983). 

10.  Techniques  to  improve  hydropower  releases  may  be  grouped  into 
three  general  areas:  forebay,  tailwater,  and  in-structure  systems. 
Examples  of  techniques  in  each  area  are: 

a.  Forebay  systems. 

(1)  Hypolimnetic  aeration/oxygenation. 

(2)  Artificial  destratification. 

(3)  Localized  mixing/local  destratification. 

(4)  Selective  withdrawal. 

b.  Tailwater  systems. 

(1)  Diffused  air  aeration. 

(2)  Weirs  or  channel  steps. 

(3)  Surface  aeration. 

(4)  Molecular  oxygen  injection. 

(5)  Miscellaneous  methods. 

c •  In-structure  systems . 

( 1 )  Air  .spiration. 

(2)  Air  or  oxygen  injection. 

Techniques  in  each  of  these  areas  will  be  discussed  in  the  following 
sections. 


Forebay  Systems 


Hypolimnetic  aeration/oxygenation 

11.  Many  forebay  systems  improve  release  DO  concentrations  by 
simply  increasing  the  DO  concentration  of  water  in  the  forebay  area. 

One  such  method,  commonly  termed  hypolimnetic  aeration  (or  oxygenation), 
involves  injecting  air  into  the  hypolimnion.  The  purpose  of  hypo¬ 
limnetic  aeration  is  to  increase  the  DO  concentrations  of  hypolimnetic 
waters  while  maintaining  the  existing  thermal  stratification. 

12.  Three  major  categories  of  hypolimnetic  aeration  have  been  sug¬ 
gested  by  Fast  and  Lorenzen  (1976):  mechanical  aeration,  air  injection, 
and  oxygen  injection.  Mechanical  aeration  has  proven  to  be  the  most 
efficient  means  of  hypolimnetic  aeration  for  shallow  lakes.  This  method 
transports  water  to  the  surface  where  it  is  mechanically  agitated  and 
returned  to  the  hypolimnion.  Air  injection  systems  mix  hypolimnetic  and 
epilimnetic  water  with  an  aerator  and  then  return  the  enhanced  waters  to 
the  hypolimnion.  Air  injection  aerators  have  demonstrated  the  highest 
oxygen  transfer  efficiencies  per  unit  energy  expended  in  pumping. 

Oxygen  injection  systems  represent  the  third  class  of  hypolimnetic 
aeration.  In  this  operation,  molecular  oxygen  is  injected  into  the 
hypolimnion  instead  of  air.  The  design  and  operation  of  hypolimnetic 
oxygenation  systems  are  discussed  in  Holland  and  Tate  (1984). 

13.  The  benefits  of  maintaining  the  thermal  characteristics  of  a 
reservoir  while  aerating  hypolimnetic  waters  include:  (a)  increasing 
the  pH  of  hypolimnetic  water  by  lowering  concentrations  of  iron,  manga¬ 
nese,  and  hydrogen  sulfide;  (b)  preventing  anoxic  conditions  which  are 
potentially  hazardous  to  fish;  (c)  maintaining  coldwater  (hypolimnetic) 
resources;  and  (d)  maintaining  heterogeneous  resources  of  water  quality 
within  the  reservoir  which  can  be  used  by  selective  withdrawal  tech¬ 
niques  to  meet  downstream  water  quality  objectives. 

14.  When  used  in  conjunction  with  hydropower  projects,  hypolim¬ 
netic  aeration  does  not  decrease  the  efficiency  of  turbines  (Speece 

et  al.  1977,  Merritt  and  Leggitt  1981).  A  major  drawback  of  air/oxygen 
injection  in  the  hypolimnion  is  that  the  volume  of  water  requiring 


enhancement,  in  most  cases,  necessitates  a  large  aeration/oxygenation 
system  (Speece  1975a).  If  oxygen,  rather  than  air,  is  required  for 
injection,  the  costs  of  oxygen  purchase,  transport,  and  storage  may  be 
quite  high  (Speece  1975b).  Hypolimnetic  aeration  may  also  result  in 
nitrogen  supersaturation,  which  can  be  hazardous  to  fish  surfacing 
downstream  of  the  project. 

Artificial  destratification 

15.  Artificial  destratification  is  another  alternative  that 
can  be  used  to  enhance  in-reservoir  and  release  water  quality  (Dortch 
1979,  Dortch  and  Holland  1980,  Fast  and  Hulquist  1982,  Holland  and 
Dortch  1984).  Destratification  requires  the  addition  of  sufficient 
energy  to  a  reservoir  to  overcome  the  buoyant  forces  associated  with 
density  stratification  and  thereby  remove  the  inhibition  to  reservoir 
circulation.  Total  reservoir  circulation  enables  the  transport  of 
oxygen  from  the  atmosphere  throughout  the  reservoir  by  convection  and 
diffusion. 

16.  Two  methods  of  artificially  destratifying  a  lake  or  reser¬ 
voir  are:  mechanically  pumping  water  (hydraulic)  and  bubbling  air 
(pneumatic).  The  hydraulic  method  jets  water  from  one  region  (hypo- 
limnion)  of  the  reservoir  into  another  region  (epilimnion) .  Pneumatic 
destratification  results  from  mixing  caused  by  an  air/water  plume  as  it 
rises  from  the  bottom  to  the  surface.  Although  pneumatic  destratifica¬ 
tion  has  been  more  widely  applied,  laboratory  tests  have  indicated  that 
hydraulic  destratification  is  possibly  more  efficient  than  pneumatic 
mixing  (Dortch  1979). 

17.  A  detrimental  effect  of  destratification,  however,  is  that  it 
results  in  the  loss  of  coldwater  resources  that  may  be  of  concern  to 
reservoir  fisheries  and  a  project's  ability  to  meet  release  coldwater 
temperature  objectives.  The  redistribution  of  thermal  energy  may  have 
the  effect  of  increasing  the  total  heat  content  of  the  reservoir.  The 
redistribution  of  nutrients  common  to  the  hypolimnion  throughout  the 
water  column  may  result  in  changes  to  the  biological  and  chemical  prop¬ 
erties  of  a  reservoir.  More  details  are  provided  in  Pastorok,  Lorenzen, 
and  Ginn  (1982). 
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Localized  mixing/ 
local  destratification 

18.  Localized  mixing  systems  are  designed  to  destratify  the  reser¬ 
voir  in  the  vicinity  of  the  outlet  as  opposed  to  total  reservoir  destra¬ 
tification  (Garton  and  Rice  1974,  Garton  and  Jarrell  1976,  Dortch  and 
Wilhelms  1978,  and  Holland  1984).  A  downward  vertical  jet  of  epi- 
limnetic  water  transports  better  quality  water  into  the  withdrawal  zone 
of  the  outlet  in  the  hypolimnion.  A  portion  of  the  transported  epilim- 
netic  water  will  then  be  withdrawn  from  the  reservoir  along  with  a  quan¬ 
tity  of  hypolimnetic  water,  thus  diluting  the  hypolimnetic  outflow  and 
improving  release  quality.  The  jet  of  water  from  the  epilimnion  may  be 
generated  by  a  number  of  techniques,  ranging  from  an  axial  flow  pro¬ 
peller  to  a  surface  pump.  This  technique  is  generally  suited  for 
smaller  flow  rates  since  a  practical  limit  of  epilimnetic  pumping  exists 
which  would  generally  provide  little  quality  enhancement  for  large- 
volume  hydropower  releases.  One  drawback  of  localized  mixing  is  acci¬ 
dental  total  lake  destratification.  A  second  is  the  likely  warming  of 
release  waters  due  to  the  increased  epilimnetic  contribution  to  the 
flow. 

19.  Local  destratification  using  a  rising  bubble  plume  has  been 
employed  as  an  alternative  to  localized  mixing  described  in  the  previous 
paragraph  (US  Army  Engineer  District  (USAED),  Savannah  1969;  Tennessee 
Valley  Authority  (TVA)  1984).  The  objective  of  this  system  is  to  create 
sufficient  mixing  locally  such  that  the  water  in  the  immediate  vicinity 
of  the  outlet  is  destratified,  resulting  in  improved  release  water 
quality.  The  disadvantages  of  the  local  destratification  system  are 
similar  to  those  for  localized  mixing. 

Selective  withdrawal 

20.  The  addition  of  selective  withdrawal  capabilities  to  a  hydro- 
power  project  is  another  alternative  that  can  potentially  improve  the 
quality  of  release  water.  Selective  withdrawal  (Wilhelms  1985)  imple¬ 
ments  the  concept  of  withdrawing  water  from  different  levels  in  a  strat¬ 
ified  reservoir  to  achieve  a  desired  release  characteristic.  For 
example,  if  warm  releases  are  desired,  surface  withdrawal  would  be  in 


order.  If  cold  water  is  desired  for  release,  a  low-level  outlet  would 
be  operated  to  withdraw  bottom  water. 

21.  A  separate  structure  may  be  added  to  the  project  (Maynord  and 
Tate  1933),  or  facilities  (or  structural  modifications)  may  be  added  to 
the  existing  structure  (George,  Dortch,  and  Tate  1980)  to  provide  con¬ 
trol  of  or  modify  the  outlet  elevation  through  which  release  water  is 
withdrawn.  Multilevel  withdrawal  provides  the  capability  to  release 
water  from  two  elevations  in  the  pool.  Thus,  control  of  release  quality 
can  be  achieved.  It  must  be  noted,  however,  that  conflicting  objectives 
may  occur  if  cold  water  with  high  DO  is  desired.  Usually,  only  surface 
water  is  high  in  DO.  Therefore,  withdrawal  of  surface  water  to  improve 
DO  would  also  increase  the  temperature  of  release  water,  resulting  in  a 
potential  conflict  with  a  coldwater  temperature  objective. 

Tailwater  Systems 


Diffused-air  aeration 

22.  One  method  of  improving  the  DO  of  the  water  in  the  tailrace  is 
by  diffused-air  aeration.  In  this  technique,  air  is  injected  into  the 
flowing  water  via  horizontal  pipes,  hoses,  or  mats  in  such  a  manner  that 
significant  air-to-water  contact  is  achieved.  This  allows  the  greatest 
efficiency  of  oxygen  absorption  from  rising  air  bubbles  into  the  sur¬ 
rounding  water.  The  effectiveness  of  such  systems  depends  upon  the 
ability  to  inject  sufficient  air  relative  to  the  water  flow  rate. 

Weirs  or  channel  steps 

23.  Significant  oxygen  uptake  may  occur  in  river  reaches  con¬ 
taining  overfall  weirs  or  multistage  channel  steps.  Oxygen  transfer 
occurs  both  from  molecular  diffusion  and  turbulent  mass  transfer.  The 
major  factor  in  determining  the  amount  of  reaeration  at  these  river 
features  is  the  fall  head  (Gameson  1957).  Additionally,  the  efficiency 
of  reaeration  through  these  structures  is  highly  dependent  on  the  depth 
of  flow.  Thus,  for  the  very  high  flows  which  are  characteristic  of  most 
hydropower  projects,  the  oxygen  transfer  would  be  limited.  Further, 
since  the  height  of  fall  determines  the  amount  of  reaeration,  the 
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addition  of  this  type  structure  to  an  existing  project  would  cause  a 
significant  loss  in  the  head  available  for  power  production. 

Surface  aeration  j 

24.  Surface  aerators  have  been  used  only  to  a  limited  extent  for 
improving  tailwater  quality.  This  alternative  usually  involves  spraying 
a  portion  of  the  water  up  into  the  air  while  agitating  and  mixing  the 
surface  waters.  The  efficiency  of  this  method  depends  upon  the  degree 
of  atmospheric  exposure  of  the  body  of  water  and  increases  with  greater 
surface  flow  velocities.  Surface  aerators  suffer  from  the  following 
drawbacks:  (a)  aesthetically  unattractive  and  noisy,  (b)  navigation 
restrictions,  and  (c)  limited  effectiveness  for  large  volumes  or  flows. 

Molecular  oxygen  injection 

25.  The  use  of  molecular  oxygen  has  been  proposed  as  an  alter¬ 
native  method  of  tailwater  aeration.  Oxygen  concentration  in  oxygen 
gas  is  about  five  times  that  of  air,  which  leads  to  a  higher  uptake 
efficiency.  It  has  been  proposed  that  fine  nozzles  or  diffuser  mats  be 
used  to  generate  the  oxygen  bubble  plume.  The  microbubbles  enhance 
oxygen  absorption  because  of  the  large  water/bubble  inter facial  area. 

The  major  problem  with  using  molecular  oxygen  is  the  cost  of  manufactur¬ 
ing  or  purchasing,  transport,  and  storage. 

Miscellaneous  methods 

26.  Miscellaneous  aeration  methods  include  Venturi  nozzles,  shaft 
aerators,  spray  cones,  and  various  surface  and  submerged  aerators  (Bohac 
et  al.  1983).  Most  of  these  methods  were  developed  for  use  at  cooling- 
water  and  sewage-effluent  outlets.  These  situations  have  small  flow 
rates  compared  to  most  hydropower  facilities  and  thus  have  limited 
applicability  downstream  of  most  CE  projects. 

In-Structure  Systems 


Air  aspiration  systems 

27.  Aspiration  systems  take  advantage  of  the  hydrodynamic/ 
hydraulic  properties  of  the  turbine  which  create  low-pressure  regions 
downstream  of  the  turbine  blades.  These  subatmospheric  pressures,  when 


vented  to  the  atmosphere,  cause  air  to  be  drawn  into  the  water  flow. 

The  existence  and  magnitude  of  a  subatmospheric  pressure  in  the  draft 
tube  are  dependent  upon  the  operating  conditions,  the  flow  rate, 
geometric  properties  of  the  turbine  and  draft  tube,  and  headwater  and 
tailwater  elevations.  The  air  flow  rate  into  the  water  flow  is  a 
function  of  the  pressure  differential  between  the  atmosphere  and  the 
draft  tube  and  the  losses  in  the  aeration  supply  line. 

28.  Francis  turbines,  the  type  used  at  many  CE  projects,  are 
usually  vented  to  the  atmosphere  (Figure  1)  during  low-flow  operations 
to  alleviate  negative  pressures  that  promote  cavitation.  Venting  under 
these  operating  conditions  also  makes  the  turbine  run  more  efficiently. 
The  automatic  venting  system  (cam-operated  valve)  that  allows  air  to  be 
drawn  into  the  draft  tube  is  called  the  "vacuum-breaker"  system. 

Usually,  the  vacuum-breaker  system  operates  over  the  lower  half  of  the 
range  of  turbine  discharges.  The  cam  action  closes  the  vacuum-breaker 
valve  at  about  a  50-percent  opening  on  the  control  gates.  Even  by 
overriding  the  automatic  closure  of  the  vacuum-breaker  valve  and  holding 
the  valve  open,  air  aspiration  is  not  extended  significantly  into  the 
higher  turbine  discharges  because  negative  gage  (subatmospheric) 
pressures  in  the  draft  tube  may  not  exist  at  higher  discharges. 


MIXTURE 

Figure  1.  Vacuum-breaker  venting  system 

29.  The  vacuum-breaker  systems  were  not  designed  to  transport 
large  volumes  of  air.  Usually,  the  piping  system  that  provides  air  to 
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the  turbine  is  long,  with  bends,  elbows,  and  valves  which  cause  a 
significant  loss  of  energy  as  air  flows  through  the  venting  network. 
This  results  in  limitations  on  the  flow  rate  of  air  that  can  be  vented 
into  the  release  flow.  Modifications  to  existing  venting  system  air- 
supply  lines  have  increased  air  flow  rates  and  enhanced  gas  transfer. 
One  example  of  a  modified  supply  line  consists  of  a  smooth  bell-mouth 
intake  that  bypasses  the  vacuum-breaker  system,  as  shown  in  Figure  2. 


AIR  WATER 
MIXTURE 

Figure  2.  Large-diameter  bell-mouth  air  intake 

30.  Generally  at  higher  flow  rates,  the  absence  of  low  pressure 
in  the  draft  tube  prevents  the  natural  aspiration  of  air.  The  deflector 
plate  concept  (Raney  1973)  was  developed  to  create  or  enhance  negative 
gage  pressures  at  the  aeration  ports  on  the  turbine  hub.  Deflector 
plates  are  placed  in  the  water  flow  upstream  of  port  openings  on  the 
turbine  hub  in  such  a  manner  that  flow  separates  from  the  turbine  hub 
over  the  aeration  portal.  This  separation  creates  low  pressures  at  the 
aeration  portal.  Even  when  free  stream  gage  pressure  is  positive,  the 
deflector  plates  create  negative  pressure  at  the  aeration  port.  The 
deflector  plate  also  increases  the  magnitude  of  the  negative  gage  pres¬ 
sures  at  the  aeration  port  when  negative  free  stream  pressures  occur 
during  lower  flow  rates.  The  larger  the  negative  gage  pressure  created 
by  the  deflector,  the  greater  the  rate  of  aspiration.  Typical  deflector 
design  and  location  on  the  hub  of  a  Francis  turbine  are  shown  in 
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a.  Daftector  plat*  dasign 


c.  Distribution  of  bafflot  on  hub 


Figure  3-  Typical  deflector  design  and  placement 


31.  The  installation  of  deflector  plates  significantly  enhances 
air  aspiration,  but  increased  hydraulic  losses  are  incurred.  The  conse¬ 
quences  of  these  losses  are  small  reductions  in  turbine  efficiency  (at 
all  operating  levels)  and  a  reduction  in  output  capacity  of  the  tur¬ 
bine.  These  losses  translate  to  lost  revenue.  Thus,  in  many  cases, 
turbine  efficiency  and  DO  enhancement  are  competing  interests.  Ideally, 
a  deflector  would  be  designed  and  placed  to  result  in  a  minimum  hydrau¬ 
lic  head  loss  while  maximizing  the  rate  of  aspiration.  Because  of  the 
additional  head  loss  caused  by  deflector  plates,  it  could  be  advan¬ 
tageous  from  a  power  generation  standpoint  to  remove  the  deflector 
plates  during  periods  when  DO  enhancement  is  not  needed. 

32.  Aspiration  can  also  be  induced  downstream  of  the  turbine  in 
the  draft  tube.  A  manifold  ring  such  as  that  shown  in  Figure  4, 
attached  to  the  periphery  of  the  draft  tube  liner,  has  been  used  to 
create  or  enhance  negative  pressures  in  the  draft  tube  (TVA  1982)  for 
essentially  the  entire  range  of  turbine  operation.  Vent  holes  on  the 
downstream  side  of  the  ring  allow  the  aspiration  of  air  into  the  release 
flow.  Uniform  spacing  and  appropriate  sizing  of  the  verit  holes  permit 
the  uniform  distribution  of  air  around  the  ring. 

33.  The  general  advantage  of  all  in-structure  aeration  techniques 
compared  to  forebay  or  tailwater  systems  is  that  all  the  outflow  from  a 
project  must  pass  through  a  confined  region.  If  the  water  quality  can 
be  enhanced  in  this  area,  it  will  impact  the  entire  outflow  from  the 
project . 

Air/oxygen  injection 

34.  The  second  method  of  in-structure  aeration  uses  an  outside 
power  source,  such  as  compressor  systems,  to  overcome  the  naturally 
occurring  hydrodynamic  pressures  in  the  turbine  to  aerate  the  hydropower 
release  water.  In  some  instances,  forced-air  injection  may  be  more 
attractive  than  induced  aspiration  since  hydraulic  losses  due  to 
deflectors  are  not  experienced.  Forced-air  venting  systems  have  the 
potential  to  aerate  flows  under  all  operating  conditions.  The  rate  of 
aeration  can  be  varied  to  correspond  with  the  degree  of  enhancement 


desired.  Compressor  systems  can  also  be  installed  to  inject  air  into 


Figure  4.  Schematic  of  ring  baffles 
installed  in  draft  tube 

the  penstock,  resulting  in  longer  contact  times  between  the  air  bubbles 
and  release  water,  thereby  improving  DO  uptake.  The  disadvantages  of  a 
forced-air  system  are:  (a)  high  initial  cost  of  equipment  (high 
volume/high-pressure  compressors  are  quite  expensive),  and  (b)  the 
operation  and  maintenance  costs  of  the  compressor  and  delivery  system. 

A  schematic  of  a  general  compressor  system  is  shown  in  Figure  5. 

Recommended  Technique 

35.  As  indicated  in  previous  paragraphs,  Bohac  et  al.  (1983)  per¬ 
formed  an  extensive  literature  review  on  the  alternatives  available  for 
improving  the  quality  of  water  released  from  hydropower  projects.  More 
particularly,  they  identified  many  efforts  to  aerate  release  flow  to 


improve  DO.  It  is  not  the  purpose  of  this  report  to  detail  their 
findings,  but  rather  to  use  their  information  to  more  clearly  define  the 
techniques  that  are  applicable  to  CE  projects.  Refinement  of  available 
information  and  data  about  a  particularly  attractive  technique  into 
guidance  for  evaluation  and  engineering  design  would  then  be  in  order. 

36.  No  single  hydropower  aeration  system  is  universally  pre¬ 
ferred.  In  general,  for  CE  impoundments,  forebay  or  in-structure 
aeration/oxygenation  systems  would  be  superior  to  tailwater  systems. 
Because  of  the  usually  large  discharges  from  a  CE  hydropower  project, 
tailwater  systems  have  limited  applicability.  In  addition  to  improving 
release  water  quality,  forebay  systems  usually  improve  in-lake  water 
quality. 

37.  The  Bureau  of  Reclamation  routinely  employs  pneumatic  de¬ 
stratification  to  maintain  in-lake  water  quality  in  water  supply  reser¬ 
voirs.  The  CE  has  employed  a  localized  pneumatic  destratification  at 
Lake  Allatoona,  Georgia  (USAED,  Savannah  1969)  for  several  years.  A 
hypolimnion  oxygenation  system  has  been  installed  at  Lake  Richard  B. 
Russell,  Georgia/South  Carolina,  and  commenced  operation  in  spring  1985 
(USAED,  Savannah  1981a,  1981b,  1982).  However,  these  systems  may 
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require  substantial  size  and  capital  outlay  because  of  the  large  volume 
of  in-lake  water  that  must  be  improved. 

38.  Based  on  the  review  of  Bohac  et  al.  (1983)  and  prototype 
studies  by  Raney  and  Arnold  (1973)  and  Mauldin  (1982),  turbine  venting 
appears  to  be  one  of  the  most  generally  applicable  techniques.  Most  CE 
projects  are  relatively  high  head  and  employ  Francis-type  turbines. 

These  turbines  (as  well  as  the  propeller  type)  are  suitable  for 
venting.  Thus,  this  technique  could  potentially  be  applied  at  many 
projects. 

39.  The  TVA  and  the  Alabama  Power  Company  (APC)  have  tested  tur¬ 
bine  venting  at  several  hydropower  projects.  Deflector  plates  for  aspi¬ 
rating  air  and  blowers  or  compressors  for  forced-air  injection  were 
tested  with  varying  degrees  of  success.  Aspiration  through  the  vacuum- 
breaker  system  normally  resulted  in  relatively  low  DO  enhancement.  In 
most  instances,  increases  in  DO  were  less  than  2  mg/d,  even  with  the 
vacuum  breaker  system  blocked  open.  Deflector  plates  increased  the  air 
flow  rate  into  the  water  discharge  and  thereby  increased  the  oxygen 
absorption.  When  the  vacuum-breaker  system  was  bypassed,  oxygen  uptakes 
of  2.5  to  3.5  mg/d  with  deflectors  were  observed.  With  a  draft  tube 
manifold  ring,  as  discussed  in  paragraph  32,  DO  improvements  of  3-5  mg/d 
were  reported. 

40.  As  discussed  earlier,  a  loss  in  turbine  efficiency  generally 
occurs  when  large  volumes  of  air  are  injected  or  aspirated  into  the 
water  flow.  At  Norris  Dam,  tests  by  the  TVA  indicated  a  loss  of  about 
3  percent  in  turbine  efficiency  when  deflector  plates  were  used  to 
aspirate  air  at  a  flow  rate  equal  to  3  percent  of  the  water  flow  rate. 

At  full-gate  operation,  output  capacity  dropped  by  about  5.5  percent. 
Similar  observations  were  made  at  test  installations  at  APC  hydropower 
facilities.  Efficiency  losses  on  the  order  of  2  percent  were  observed 
in  conjunction  with  the  installation  of  an  aeration  manifold  ring. 

These  losses  are  essentially  a  head  loss;  consequently,  the  percent  loss 
associated  with  these  alternatives  may  tend  to  decrease  with  increasing 
head. 

41.  The  synthesis  of  design  procedures  to  implement  the  turbine 


venting  alternatives  (either  aspiration  or  injection)  is  a  necessary 
step  in  the  application  of  turbine  venting  technology.  Therefore,  the 
physical  processes  that  impact  reaeration  during  turbine  venting  must  be 
identified.  Descriptions  of  these  processes  must  be  developed  so  that 
an  accurate  prediction  of  oxygen  uptake  can  be  made.  Thus,  an  evalu¬ 
ation  of  the  effectiveness  of  turbine  venting  can  be  accomplished  for  a 
specific  stratification,  discharge,  and  outlet  geometry  situation. 

42.  The  relationship  of  turbine  performance  and  turbine  venting 
must  be  understood  to  determine  the  economic  impacts  of  turbine 
venting.  The  remainder  of  this  report  presents  the  effort  to  further 
understand  the  processes  and  relationships  that  govern  oxygen  uptake 
with  turbine  venting  and  develop  descriptions  of  those  that  can  be  used 
for  evaluation  and  design  guidance. 


PART  III:  FIELD  STUDIES  FOR  PROCESS  DESCRIPTION 


Field  Study  Site 

43.  As  stated  in  previous  paragraphs,  adequate  design  guidance  to 
implement  a  turbine  venting  system  has  not  been  completely  developed. 
Hence,  field  studies  were  initiated  to  quantify  the  design  aspects  and 
more  clearly  define  the  impacts  of  turbine  venting.  Field  studies  were 
conducted  at  the  Clarks  Hill  Dam  on  the  Savannah  River.  The  project  has 
seven  Francis-type  turbines  rated  at  a  power  capacity  of  40  MW  each. 
Clarks  Hill  Dam  is  a  peaking  power  project;  thus,  the  turbines  are  usu¬ 
ally  operated  daily  for  6  to  10  hr  at  levels  from  35-  to  100-percent 
capacity. 

44.  The  vacuum-breaker  venting  system  (described  in  paragraph  28) 
at  this  site  consisted  of  a  cam-operated  valve  and  an  8-in.*  air-supply 
pipe.  The  vacuum-breaker  system  operated  at  wicket  gate  openings  below 
43  percent.  On  the  turbine  hub,  there  were  eight  equally  spaced  3-  by 
7-in.  holes  through  which  the  vented  air  entered  the  flow. 

45.  Two  turbine  venting  techniques  were  investigated :  (a)  aspi¬ 

ration  and  (b)  injection.  Aspiration  was  investigated  on  Unit  2  (no 
deflector  plates)  with  the  existing  vacuum-breaker  system  forced  open 
for  the  entire  range  of  wicket  gate  settings.  A  second  aspiration 
arrangement  was  tested  by  replacing  the  vacuum-breaker  valve  with  a 
10-in. -diam  smooth  bell-mouth  intake.  A  final  series  of  tests  was 
conducted  on  Unit  2  in  which  air  was  forced  into  the  turbine  with  a  low- 
pressure  (approximately  0.5  psi,  although  a  high-pressure  compressor 
would  have  proven  more  effective)  25-hp  blower. 

46.  Unit  4  was  fitted  with  deflector  plates  to  decrease  the  pres¬ 
sure  at  the  vent  holes  on  the  turbine  hub.  The  deflector  plate  design 
was  based  on  prior  work  by  TVA  (1981),  Mauldin  (1982),  and  the  APC 
(Raney  1973,  1975).  The  deflectors  on  Unit  4  were  constructed  of 
stainless  steel  with  a  45-deg  leading  angle.  The  deflector  plate  was 

*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 


8  in.  wide,  3  in.  deep,  and  3  in.  high  (Figure  6).  The  deflectors  were 
oriented  and  welded  in  place  on  the  hub  about  4  in.  upstream  of  the  vent 
holes  so  that  the  lowest  pressure  in  the  wake  of  the  deflector  occurred 
over  the  vent  holes  (Mauldin  1982).  Tests  were  conducted  with  the 
vacuum-breaker  system  blocked  open.  The  second  configuration  tested  on 
Unit  4  was  the  replacement  of  the  vacuum-breaker  valve  with  a  bell-mouth 
intake  to  increase  the  air  flow  into  the  turbine. 


47.  Both  units  were  tested  under  "no-air"  conditions,  i.e.,  vent¬ 
ing  mechanisms  held  shut,  to  establish  baseline  performance  data  on  tur¬ 
bine  operation  and  oxygen  absorption.  The  operating  threshold  at  Clarks 
Hill  powerhouse  was  30  percent.  Power  production  limitations  prevented 
gate  settings  greater  than  85  percent.  Therefore,  wicket  gate  settings 
ranged  from  30  to  85  percent  on  all  turbine  venting  tests. 

Data  Collection  and  Instrumentation 

48.  To  develop  engineering  guidance  for  turbine  venting,  the 
processes  that  affect  oxygen  uptake  and  turbine  performance  must  be 
understood.  Further,  for  rigorous  analyses,  mathematical  descriptions 


of  these  processes  have  to  be  developed.  By  examining  the  effects  of  a 
process,  often  the  cause  can  be  identified.  Thus,  field  data  collection 
and  analysis  can  produce,  when  coupled  with  appropriate  theoretical 
analyses,  the  mathematical  relationships  for  the  processes  that  cause 
the  effects.  A  greater  understanding  of  those  processes  is  thereby 
acquired . 

49.  The  enhancement  of  DO  concentrations  in  hydropower  releases 
is  the  objective  of  turbine  venting.  Therefore,  upstream  (penstock)  and 
downstream  (tailrace)  DO  concentrations  were  measured. 

50.  When  gas  transfer  occurs  due  to  air  bubbles  in  the  water 
body,  more  air  bubbles  should  result  in  more  gas  transfer.  Thus,  for 
turbine  venting,  the  flow  rate  of  air  into  the  turbine  significantly 
impacts  the  reaeration  (oxygen  uptake)  process.  The  air  flow  into  the 
turbine  was  determined  for  all  the  conditions  tested.  Additionally, 
since  the  air  flow  rate  is  affected  by  the  pressure  in  the  turbine  and 
draft  tube,  the  turbine  head  cover  pressure  was  measured. 

51.  For  any  hydraulic  situation  when  air  is  entrained  or  injected 
into  a  pressure  flow  (draft  tube)  condition,  the  potential  exists  for 
gas  supersaturation  to  occur.  Total  dissolved  gas  pressure,  from  which 
the  level  of  saturation  can  be  determined,  was  also  measured.  For  com¬ 
plete  analysis  of  dissolved  gas  saturation  and  to  determine  the  level  of 
nitrogen  saturation,  DO,  water  temperature,  and  barometric  pressure  data 
are  required. 

52.  Power  output,  water  flow  rate,  and  headwater  and  tailwater 
elevations  were  collected  to  evaluate  the  impact  of  injected  or  aspi¬ 
rated  air  on  turbine  operation.  These  turbine  performance  data  are  re¬ 
quired  to  calculate  the  operating  efficiency  of  the  turbine.  Wicket 
gate  settings  which  indicate  the  level  of  turbine  operation  were  re¬ 
corded.  From  these  data,  a  comparison  of  the  various  aeration  tech¬ 
niques  and  associated  impacts  on  turbine  operation  can  be  performed. 

53.  Data  were  collected  from  several  locations  at  the  dam.  The 
tailrace  (downstream)  data  collection  point  was  located  in  the  turbine 
release  downstream  of  the  aeration  bubble  plume.  Simultaneous  data  col¬ 
lection  at  each  of  the  collection  points  was  accomplished  by  use  of 


portable  radio  communication.  Table  1  lists  the  data  types  collected  in 
these  field  studies  and  describes  the  measurement  techniques  used  in 
their  collection. 


Table  1 

Data  and  Method  of  Calculation 


_ Parameter _ 

Penstock  and  DO  (mg/t) 

Water  temperature  (°C) 
Barometric  pressure  (mb) 
Air  flow  rate  (cfs) 


Head  cover  pressure  (ft  water) 

Total  dissolved  gas  pressure 
(TDGP)  (mm  HG) 

Power  output  (MW) 

Water  flow  rate  (cf3) 


Headwater  and  tailwater 
elevations  (ft) 

Wicket  gate  opening  (%  open) 


_ Method  of  Collection _ 

Yellow  Springs  Instrument  Co.  (YSI)  polarographic  DO 
probes  (calibrated  with  the  Modified  Winkler 
method) . 

Water  temperature  was  measured  with  mercury  thermom¬ 
eters  and  the  YSI  probe. 

Readings  were  made  three  times  daily  with  a  portable 
barometer. 

An  air-velocity  profile  was  measured  with  a  hot-wire 
anemometer,  or  a  Pitotstatic  tube  and  manometer. 

Air  flow  rate  was  determined  by  integrating  the 
point  velocity  readings  over  the  area  of  the 
condu i t . 

A  pressure  gage  was  read  at  the  turbine. 

A  saturometer  designed  at  WES  (Wilhelms  1984)  was 
used  to  measure  TDGP.  Total  nitrogen  concentration 
was  calculated  from  TDGP  readings  and  DO  readings  by 
assuming  all  gas  in  the  water  that  was  not  oxygen 
was  nitrogen. 

Revolutions  of  the  watt-hour  meter  in  the  control 
room  were  counted  over  timed  intervals. 

Readings  were  taken  from  a  differential  manometer 
connected  to  Winter-Kennedy  pressure  taps  on  the 
turbine  scroll  case  and  from  a  digital  flow  rate 
meter  in  the  control  room. 

Readings  were  taken  from  meters  in  the  control  room. 


Readings  were  made  from  a  meter  in  the  control  room 
and  from  Indicators  on  the  turbine. 


Note:  Tables  A1-A7  of  Appendix  A  present  data  from  the  1982  field  studies.  Tables  A8-A9 
present  data  from  the  1981  field  studies. 


PART  IV:  RESULTS  OF  FIELD  TESTING 


Data  Reduction 


54.  To  permit  comparative  analyses  of  the  Clarks  Hill  data,  the 
power  output  and  discharge  measurements  were  adjusted  to  a  gross  head 
of  146  ft.  The  following  relationships  (Pfau  1948)  were  used  in  making 
these  adjustments. 


where 

PADJ*^ADJ  =  power  output  and  flow  rate,  respectively, 
adjusted  to  146-ft  head 

PACT'^ACT  =  observed  power  output  and  flow  rate, 
respectively 

^ADJ  =  adjusted  gross  head,  146  ft 

U 

ACT  =  observed  gross  head,  ft 

55.  For  each  condition  tested,  the  adjusted  data  were  plotted 
against  wicket  gate  opening  to  facilitate  analysis.  Figures  7a  and  7b 
show  the  power  and  discharge  graphs,  respectively,  for  Unit  2,  without 
air  injection  or  aspiration.  These  plotted  data  were  smoothed  by  hand, 
as  indicated  by  the  solid  lines.  Subsequently,  "smoothed"  power,  dis¬ 
charge,  and  wicket  gate  data  were  digitized  from  these  curves  and  were 
used  in  process  analyses.  Plots  similar  to  Figures  7a  and  7b  for  all 
the  venting  techniques  (Units  2  and  4)  are  presented  in  Figures  B1-B14 
of  Appendix  B. 

56.  Turbine  efficiency  was  calculated  using  the  following 
relationship : 


E 


(3) 


where 


KP 
=  QH 


E  =  turbine  efficiency,  percent 
K  =  conversion  factor,  ft^/sec  -  MW 
P  =  power  output,  MW 
Q  =  turbine  discharge,  cfs 

H  =  headwater  elevation-tailwater  elevation,  ft 


The  K  terra  has  a  value  of  1.18  (10^)  and  accounts  for  the  conversion 
of  units  and  the  effects  of  power  transformers.  In  this  computation, 
the  P  and  Q  values  were  digitized  from  the  smoothed  plots  of  power 
and  discharge  versus  wicket  gate  setting,  respectively. 


Analysis  of  Turbine  Performance 

57.  As  stated  previously,  air  flow  rate  was  considered  to  be  an 
important  parameter  affecting  DO  uptake.  Air  flow  rate  was  also  found 
to  significantly  impact  turbine  performance  (Raney  1973;  Buck,  Miller, 
and  Sheppard  1980;  Mauldin  1982).  Thus,  to  evaluate  that  impact,  the 
no-air  test  results  were  used  as  the  base  condition  for  analysis  of 
turbine  efficiency  changes.  Turbine  efficiency,  as  computed  with  Equa¬ 
tion  3,  was  plotted  against  wicket  gate  setting  for  the  no-air  condition 
and  for  each  of  the  venting  arrangements.  Figure  8  shows  efficiency 
versus  wicket  gate  for  the  no-air  tests  and  the  blower  tests  on  Unit  2. 
Similar  plots  were  developed  for  each  of  the  venting  techniques  on 
Units  2  and  4  and  are  shown  in  Figures  B15-B19.  This  difference  in  the 
efficiencies  (between  vented  and  no-air  conditions)  shown  on  these  plots 
is  indicative  of  the  impact  of  venting  (air  injection  or  aspiration)  on 
turbine  performance. 

58.  These  figures  indicate  that,  in  general,  air  introduction 
caused  a  small  loss  in  efficiency  compared  to  the  nonvented  condition. 
The  loss  was  particularly  pronounced  for  the  lower  (less  than  50-per¬ 
cent)  wicket  gate  openings  when  relatively  large  volumes  of  air  were 
being  drawn  or  forced  into  the  turbine.  Referring  to  Figure  8,  the  loss 
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Figure  8.  Efficiency  versus  wicket  gate,  Unit  2 

experienced  by  the  Unit  2  turbine  at  a  45-percent  wicket  gate  was  about 
4  percent.  At  a  55-percent  wicket  gate,  the  loss  was  about  1.5  percent. 
Similar  trends  are  apparent  from  the  other  venting  test  conditions  by 
examining  Figures  B15-B19. 

59.  It  should  be  noted  that  this  efficiency  loss  is  a  real  loss 
in  power  output  and  not  "power  deferred,"  i.e.,  a  decrease  in  discharge 
that  "saves"  water  for  later  power  production.  Plots  of  turbine  effi¬ 
ciency  versus  discharge  (Figures  B20-B24)  such  as  that  shown  in 

Figure  9,  clearly  demonstrate  that  the  loss  in  efficiency  was  not  offset 
by  a  corresponding  decrease  in  discharge  rate.  A  1-percent  loss  in  ef¬ 
ficiency  would  be  accompanied  by  a  1 -percent  decrease  in  discharge  rate 
if  the  "power  deferred"  concept  were  applicable  for  these  techniques. 

60.  Of  particular  interest  in  Figure  8  is  the  slight  efficiency 
increase  at  and  above  a  wicket  gate  setting  of  60  percent  for  the  vented 
condition  compared  to  the  no-air  condition.  At  these  gate  settings,  the 
air  flow  rates  were  small  relative  to  the  water  flow  rates,  compared  to 
those  at  the  lower  gate  settings  (Figure  10).  The  improved  performance 
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can  probably  be  attributed  to  smoother  turbine  operation  resulting  from 
the  introduction  of  the  small  air  volume.  In  an  analysis  of  Figures  8 
and  10,  it  can  be  concluded  that  air,  vented  into  the  turbine,  can  im¬ 
prove  or  degrade  the  performance  of  the  turbine  depending  upon  the 
amount  of  air. 

61.  To  permit  a  more  general  analysis  of  the  effects  of  injected 
or  aspirated  air  on  the  turbine  performance,  the  air  flow  parameter  was 
nondimens ionalized  by  dividing  it  with  the  water  flow  rate.  Thus,  an 
air  flow-to-water  flow  ratio  of  0.02  means  that  air  flow  was  2  percent 
of  the  water  flow. 

62.  To  obtain  these  ratios,  observed  air  flow  (Q-air)  and  ob¬ 
served  water  discharge  (Q-water)  data  were  plotted  against  wicket  gate 
opening  in  Figure  10  for  the  Unit  2  blower  tests.  Similar  graphs  for 
each  of  the  venting  techniques  are  presented  in  Appendix  B 

(Figures  B25-B29).  Flow  rate  information  at  any  given  gate  setting  was 
interpolated  or  extrapolated  from  observed  data. 

63.  By  referring  to  the  efficiency  curves  shown  in  Figures  8  and 
B15-B19  and  the  air  flow  and  water  flow  curves  shown  in  Figure  10  and 
B25-B29,  for  a  given  wicket  gate  opening  an  air  flow-to-water  flow 
(Q-air/Q-water)  ratio  and  an  efficiency  loss  can  be  determined.  Plot¬ 
ting  efficiency  loss  against  the  ratio  of  air  flow  to  water  flow  for 
Unit  2  (Figure  11)  shows  the  change  in  turbine  performance  due  to  vent¬ 
ing  air  into  the  turbine.  In  most  tests  an  efficiency  loss  was  ob¬ 
served;  however,  at  very  small  Q-air/Q-water  ratios,  a  slight  effi¬ 
ciency  improvement  was  observed  for  Unit  2.  As  stated  earlier,  this 
was  probably  the  result  of  smoother  turbine  operation  due  to  the  intro¬ 
duction  of  very  small  air  volumes. 

64.  For  Unit  2,  without  deflector  plates,  the  efficiency  loss 
was  linearly  related  to  Q-air/Q-water  except  for  the  tests  of  the 
vacuum-breaker  system.  It  appears  that  this  nonlinearity  was  due  to  the 
flow  characteristics  of  the  vacuum-breaker  system,  i.e.,  air  flow  was 
limited  (compare  air  flow  in  Figure  B25  and  the  vacuum-breaker  effi¬ 
ciency  loss  in  Figure  11).  Figure  12  shows  efficiency  loss  plotted 
against  Q-air/Q-water  for  Unit  4  (the  turbine  with  deflectors).  This 
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Figure  11.  Efficiency  loss  versus  air 
flow-to-water  flow  ratio,  Unit  2 
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Figure  12.  Efficiency  loss  versus  air 
flow-to-water  flow  ratio,  Unit  4 
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graph  indicated  that  there  was  nearly  always  an  efficiency  loss.  This 
maybe  a  consequence  of  the  relatively  large  volumes  of  air  that  were 
aspirated  into  the  turbine  over  the  entire  range  of  wicket  gate  set¬ 
tings  (Figures  B28-B29)  as  a  result  of  the  deflectors  or  a  combined 
effect  of  the  air  and  the  hydraulics  of  the  deflectors. 

Analysis  of  Oxygen  Transfer  Characteristics 

65.  To  analyze  the  DO  uptake  due  to  turbine  venting,  an  under¬ 
standing  of  the  physics  that  affect  gas  transfer  is  necessary.  The 
force  driving  the  gas  transfer  process  is  the  difference  in  partial 
pressure  of  the  oxygen  in  the  water  and  in  the  air  at  the  air/water 
interface.  If  the  partial  pressures  of  oxygen  in  the  water  and  in  the 
air  are  equal,  then  no  driving  force  exists  and  the  net  oxygen  transfer 
between  the  air  and  the  water  is  zero.  Under  this  thermodynamically 
equilibrated  condition,  the  water  is  considered  "saturated."  If  the 
partial  pressures  in  the  water  and  in  the  air  are  unequal,  then  a  force 
exists  to  cause  an  oxygen  transfer. 

66.  A  measure  of  this  force  is  the  "saturation  deficit,"  which 
can  be  quantified  by  using  the  concentration  of  the  oxygen  in  the  water 
and  the  saturation  concentration  defined  by  Henry's  Law  (Schroeder  1977) 


C  =  kp  (4) 

o 


where 


Cg  =  saturation  concentration  for  oxygen 

k  =  proportional  coefficient  for  oxygen 

p  =  partial  pressure  of  oxygen  in  atmosphere  across 
air/water  interface 

Henry's  Law  simply  states  that  water  (at  a  given  temperature)  can  con¬ 
tain  an  amount  of  DO  that  is  linearly  proportional  to  the  partial  pres¬ 
sure  of  that  gas  in  the  atmosphere  across  air/water  interface.  The 
"saturation  deficit"  is  the  difference  between  the  saturation  concen¬ 


tration  and  the  actual  concentration  of  oxygen  in  the  water  and  is 
describe  by 


67.  Reaeration  or  oxygen  uptake  is  usually  considered  to  be  a 
first-order  reaction  that  is  described  in  in  terms  of  initial  and  final 
(or  upstream  and  downstream)  saturation  deficits.  Therefore,  the  oxygen 
transfer  characteristics  of  turbine  venting  were  analayzed  on  the  basis 
of  the  relationships  between  the  upstream  and  downstream  deficits  and 
the  air-to-water  flow  ratio.  To  nondimens ionalize  the  quantities,  a 
"deficit  ratio"  was  defined  as  the  ratio  of  the  downstream  deficit  (DO 
deficit  of  water  leaving  vicinity  of  tailrace  area)  to  the  upstream 
deficit  (DO  deficit  of  water  in  penstock). 

68.  Figure  13  shows  the  relationship  between  deficit  ratio  and 
Q-air/Q-water .  Least  squares  regression  analysis  of  these  data  resulted 
in 


^4^  =  ^  =  0.67  .  0.424  e-70'91- 
s  u  u 


(6) 


where 


Cg  =  temperature-dependent  saturation  concentration  of 
oxygen 

Cd,Cu  =  downstream  and  upstream  DO  concentrations, 
respectively 

Dd,Du  =  downstream  and  upstream  DO  deficits,  respectively 
r  =  air  flow-to-water  flow  ratio 


Figure  13  clearly  indicates  that  the  oxygen  transfer  characteristics  of 
turbine  venting  are  a  function  of  Q-air/Q-water  .  The  scatter  in  the 
data  suggests  that  other  conditions  also  impact  oxygen  uptake. 

69.  Figure  13  and  Equation  6  suggest  that  there  is  a  maximum  re¬ 
duction  in  deficit  ratio  which  can  be  achieved  with  turbine  venting. 

The  exponential  term  in  the  equation  becomes  very  small  as  the  air  flow- 
to-water  flow  ratio  increases.  This  indicates  that  with  a  large  r  , 
the  downstream  deficit  would  be  about  70  percent  of  the  upstream  de¬ 
ficit.  For  example,  if  the  upstream  deficit  was  8.0  mg/e,  the 
downstream  deficit  would  be  approximately  5.6  mg/e  (Df/Di  =  0.70)  if  the 
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Figure  13.  Deficit  ratio  versus  air  flow-to-water 
flow  ratio  (Q-air/Q-water) 

air  flow  was  3.7  percent  of  the  water  flow  (Q-air/Q-water  =  0.037).  A 
DO  uptake  of  2.4  mg/l  would  be  experienced  under  these  condition. 
Further,  Equation  6  and  Figure  13  indicate  that  even  with  6-percent  air 
flow,  the  oxygen  uptake  would  increase  to  only  2.6  mg/i.  Thus,  further 
increases  in  air  flow  may  produce  only  marginal  impacts  on  the  DO 
concentration  of  the  release. 

Comparison  of  Techniques 

70.  A  comparison  of  the  venting  systems  for  each  of  the  turbines 
is  presented  in  Table  2.  Wicket  gate  settings  of  40,  50,  and  60  percent 
were  selected  as  the  settings  at  which  comparisons  could  be  made.  At 
these  gate  openings,  all  the  conditions  (no-air,  vacuum  breaker,  bell- 
mouth  intake,  and  blower)  were  tested  on  Unit  2,  with  the  exception  of 
the  blower  at  40-percent  gate  opening;  all  the  techniques  except  the 
blower  were  tested  on  Unit  4.  To  compare  the  oxygenation  capabilities 
of  each  system,  the  percent  of  deficit  satisfied  is  shown  in  Table  2. 


The  deficit  satisfied  is  the  difference  between  the  upstream  and  down¬ 
stream  deficits  divided  by  the  upstream  deficit. 


SD  = 


=  1  - 


(7) 


where  SD  =  deficit  satisfied. 

71.  Referring  to  the  example  discussed  in  paragraph  69,  if  the 
downstream  deficit  was  70  percent  of  the  upstream  deficit,  then 


and 


0.70 


SD  =  1  -  0.70  =  0.30 


(8) 


or  30-percent  deficit  satisfaction  was  achieved.  Efficiency  change  is 
also  shown  in  Table  2.  These  values  were  computed  by 


AE  =  E  -  E 
o  c 


(9) 


where 


AE  =  efficiency  change,  percent 

EQ  =  efficiency  of  turbine  under  no-air  condition,  percent 

Ec  =  Efficiency  of  turbine  under  each  venting  technique, 
percent 

Although  the  values  presented  in  Table  2  represent  individual  test 
results,  they  provide  a  basis  for  comparing  the  systems  in  general 
terms. 

72.  From  the  table,  the  trade-offs  between  efficiency  and  gas 
transfer  become  obvious.  Generally,  the  higher  efficiency  losses  are 
associated  with  larger  deficit  satisfaction  (more  oxygen  uptake).  For 
example,  on  Unit  2,  at  a  wicket  gate  setting  of  50  percent,  the  vacuum 
breaker  resulted  in  an  efficiency  loss  of  0.7  percent  with  19  percent  of 
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the  deficit  satisfied.  Replacing  the  vacuum  breaker  with  the  bell-mouth 
intake  allowed  more  air  flow,  and  the  deficit  satisfied  increased  to 
23  percent.  However,  the  efficiency  loss  increased  as  well,  to 
1 . 3  percent . 

73.  An  increase  in  turbine  efficiency  at  a  wicket  gate  setting  of 
60  percent  was  observed  on  Unit  2  for  all  systems  tested.  The  gas 
transfer  at  this  setting  was  very  small,  as  might  be  expected  for  the 
very  low  air  flow  being  vented  into  the  turbine  (see  Figures  B25-B27). 
Thus,  benefits  of  improved  efficiency  and  significant  gas  transfer  can¬ 
not  be  realized  simultaneously.  The  increased  efficiency  was  not  ob¬ 
served  for  any  test  on  Unit  4.  This  may  have  been  caused  by  the  pre¬ 
sence  of  deflectors  and/or  the  large  air  flow  rate  (Figures  B28-B29). 

The  turbine  efficiency  losses  and  gas  transfer  were  generally  greater 
with  deflectors  than  without. 

74.  Based  on  system  comparison,  use  of  the  vacuum  breaker  has  a 
moderate  impact  on  both  efficiency  and  gas  transfer.  However,  when 
coupled  with  deflectors,  the  effects  are  significantly  increased.  Using 
the  bell-mouth  intake  in  place  of  the  vacuum-breaker  system  signifi¬ 
cantly  decreased  the  efficiency  in  the  lower  range  of  wicket  gate  set¬ 
tings.  Since  these  low  wicket  gate  settings  are  generally  not  used  for 
power  generation,  these  efficiency  losses  may  be  inconsequential.  At 
higher  wicket  gate  settings,  the  moderate  increase  in  gas  transfer 
associated  with  the  bell-mouth  intake,  compared  with  the  vacuum-breaker 
system,  may  warrant  accepting  the  slight  decrease  in  efficiency.  There 
was  little  difference  between  the  effects  of  the  bell-mouth  intake  and 
blower  since  the  small  blower  (7,000  cfm  air  at  0.5  psi)  did  r.ot 
significantly  increase  the  air  flow  above  the  air  flow  rate  with  the 
bell-mouth  intake  (Figures  B26-B27). 


PART  V:  MODELING  OF  REAERATION  THROUGH  A 
VENTED  HYDROTURBINE 


75.  Development  of  a  technique  to  predict  the  oxygen  uptake  due 
to  turbine  venting  is  requisite  for  subsequent  development  of  engi¬ 
neering  design  guidance.  The  technique  must  include  in  its  formulation 
factors  that  account  for  as  many  of  the  actual  physical  processes  as 
possible.  If  the  predictive  technique  accurately  models  the  processes, 
the  potential  effects  of  a  proposed  turbine  venting  installation  could 
be  evaluated.  If  the  impact  of  turbine  venting  could  be  estimated,  the 
costs  (in  terms  of  efficiency  loss)  of  the  increased  oxygen  in  the  re¬ 
lease  could  be  determined  and  alternatives  compared.  This  application 
will  be  discussed  with  an  example  in  the  concluding  section  of  Part  V. 

76.  The  data  discussed  in  Part  IV  and  presented  in  Appendix  A 
give  us  only  the  results  of  the  various  processes  affecting  release 
DO.  Thus,  it  is  necessary  to  identify  each  process  and  the  impacts  of 
each.  Once  this  is  accomplished,  a  mathematical  description  of  each 
process  can  be  developed.  Ultimately,  the  process  descriptions  must  be 
combined  to  produce  a  complete  mathematical  model  of  the  reaeration 
effects  of  turbine  venting.  This  numerical  model  can  then  serve  as  a 
predictive  technique  with  which  proposed  turbine  venting  may  be  eval¬ 
uated.  The  following  section  shows  the  development  of  a  numerical  model 
of  reaeration  due  to  turbine  venting. 

Model  Development 

77.  Analysis  of  the  data  from  the  no-air  tests  indicated  that 
there  was  a  DO  uptake  in  the  release  in  the  tailrace  area.  When  the 
turbine  was  vented  and  air  was  introduced  into  the  flow  (by  aspiration 
or  injection),  additional  DO  uptake  was  observed.  It  was  concluded  that 
two  processes  were  causing  gas  transfer.  Turbulent  reaeration  in  the 
tailrace  was  responsible  for  the  oxygen  uptake  during  the  no-air  tests. 
The  additional  uptake,  during  venting,  was  due  to  oxygen  transfer  from 
the  injected  air  to  the  release  water.  Figure  14  shows  the  conceptual 
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Figure  14.  Conceptual  relationship  of  gas  transfer  processes 

relationship  of  these  two  processes.  For  development  of  the  predictive 
model,  the  effects  of  these  processes  were  mathematically  described  and 
superimposed  to  form  a  numerical  model  of  the  DO  uptake  due  to  turbine 
venting. 

78.  Turbulent  reaeration  in  the  tailrace  area  was  considered  a 
function  of  the  available  energy  at  the  draft  tube  exit.  The  Energy 
Dissipation  Model  (EDM)  (Tsivoglou  and  Wallace  1972,  Wilhelms  and 
Smith  1981)  was  used  to  account  for  the  DO  uptake  during  the  no-air 
tests.  The  EDM  describes  reaeration  by 
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saturation  concentrations  for  ambient  water 
temperature  T,  mg/t, 
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Cd,Cu  =  tailrace  and  penstock  DO  (without  venting),  mg/ll 
Dd,Du  =  tailrace  and  penstock  oxygen  deficits,  mg/i 
AE  =  energy  dissipation  coefficient,  ft-^ 
cT  =  energy  dissipated  in  tailrace,  ft-lb/lb 


79.  The  energy  dissipated  in  the  tailrace  can  be  approximated  by 
writing  Bernoulli's  equation  at  the  draft  tube  outlet  and  at  a  point  in 
the  tailrace  downstream  of  the  high-turbulence  area. 


(vtr) 

2g 


+  Zrpp  +  AE 


where 


PDT’PTR 


VDT’VTR  = 


ZDT,ZTR 


pressure  at  draft  tube  (DT)  center  line  and  mid¬ 
depth  in  the  tailrace  (TR),  respectively,  lb/ft^ 

specific  weight  of  water,  62.4  lb/ft^ 

average  velocity  at  DT  and  TR,  respectively,  ft/sec 

2 

acceleration  due  to  gravity,  32.2  ft/sec 

elevation  above  arbitrary  datum  for  location  at  DT 
and  TR,  respectively,  ft 


80.  The  variables  are  also  defined  in  Figure  15.  It  was  assumed 
that  the  difference  in  water  surface  elevations  between  these  two  points 
was  negligible,  thus 


It  was  further  assumed  tha  tthe  downstream  velocity  was  very  small  re¬ 
lative  to  dreft  tube  velocity,  i.e.,  ^V^R^/2g  =  0  .  Thus,  the  energy 

dissipated  in  the  tailrace  can  be  approximated  by 


(vdt) 


81.  The  energy  dissipation  coefficient  is  dependent  on  tem¬ 
perature  (Churchhill,  Elmore,  and  Buckingham  1962;  Tsivoglou  and 


Figure  15.  Variable  definitions  for  Bernoulli  equation 


Wallace  1972)  and  the  relationship  is  described  by 


cT  =  c2Q( 1.024) 


(T-20) 


where 


“T * c20  =  ener8y  dissipation  coefficient  for  water 
temperatures  of  (ambient)  and  20°  C,  ft-' 

T  =  ambient  water  temperature,  °C 


Substituting  Equations  13  and  14  into  Equation  10  results  in  the 
following  mathematical  description  of  reaeration  due  to  downstream 


turbulence. 
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82.  Using  the  observed  DO  readings  from  the  no-air  tests  and 
continuity 

„  .  Q 


'DT  "  A 


DT 


(16) 


where 


Q  =  discharge,  cfs 


*DT 


cross-sectional  flow  area  of  draft  tube  exit,  ft' 


to  determine  draft  tube  velocity,  a  least  squares  regression  analysis 
produced  a  C20  value  of  0.15  per  foot.  The  resulting  equation  can 
be  used  as  a  reaeration  model  when  venting  is  not  in  operation.  It 
must  be  remembered  that  the  turbulent  reaeration  that  occurred  in  the 
tailrace  exists  with  or  without  venting.  Thus,  when  venting  is  em¬ 
ployed,  this  oxygen  transfer  must  still  be  included,  as  shown  in  the 
conceptual  model  in  Figure  14. 

83.  As  stated  earlier,  the  release  DO  improved  significantly 
when  venting  was  initiated  compared  to  the  no-air  condition.  Thus,  the 
gas  transfer  occurring  as  a  result  of  the  vented  air  is  the  process 
that  must  be  described.  The  discussion  of  DO  uptake  (paragraphs  65-69) 
concluded  that  the  reaeration  due  to  venting  was  a  function  of  the  ratio 
of  air  flow  to  water  flow.  Ordinarily,  the  gas  transfer  process  is 
considered  to  be  a  first-order  reaction  and,  as  such,  is  mathematically 
described  by 


where 

K  =  exchange  coefficient,  sec"1 

t  =  time  of  flow  from  an  upstream  location  to  a  downstream 
location,  sec 

84.  A  linear  relationship  between  the  exchange  coefficient  and 
r  (Q-air/Q-water )  was  assumed. 


where 


a  =  coefficient  of  gas  transfer,  sec" 
r  =  ratio  of  air  flow  to  water  flow,  dimensionless 
Hence,  Equation  17  becomes 


u 


To  apply  Equation  19,  the  effects  of  hydraulic  forces  must  be  estimated 
as  well  as  the  travel  time  from  an  upstream  point  to  a  downstream  point. 

85.  The  hydraulic  forces  that  act  on  the  air  bubbles  as  they 
travel  from  the  venting  port  on  the  turbine  hub  to  the  tailrace  must  be 
understood.  Those  forces  change  the  thermodynamic  state  of  the  air 
bubbles  as  they  move  through  the  draft  tube  and  thereby  impact  the  gas 
transfer  from  the  bubbles  to  the  water.  Consider  the  thermodynamic 
state  of  an  air  bubble  at  the  surface  of  a  water  body.  The  sum  of  the 
partial  pressures  of  the  gases  that  comprise  the  bubble  is  essentially 
atmospheric  pressure;  therefore,  the  partial  pressure  of  oxygen  in  the 
bubble  is  equal  to  the  partial  pressure  of  oxygen  in  the  atmosphere. 

86.  According  to  Henry's  Law  (Equation  4),  an  oxygen  saturation 
concentration  for  the  water  surrounding  the  bubble  can  be  determined. 

If  the  oxygen  concentration  in  the  water  surrounding  the  bubble  is  equal 
to  the  saturation  concentration,  the  oxygen  gas  in  the  bubble  and  in  the 
water  is  at  thermodynamic  equilibrium  which  is  characterized  by  no  net 
oxygen  transfer  from  the  bubble  to  the  water. 

87.  If  the  air  bubble  is  forced  deeper  into  the  water  body,  the 
hydrostatic  pressure  acting  on  the  bubble  will  increase  with  an  iden¬ 
tical  increase  in  the  pressure  of  the  air  inside  the  bubble.  The  in¬ 
crease  in  bubble  pressure  also  increases  the  partial  pressures  of  the 
gases  that  make  up  the  air.  According  to  Henry's  Law,  there  would  be  a 
proportional  increase  in  the  saturation  concentration.  Ftr  example,  at 
the  surface  of  an  impoundment,  for  a  water  temperature  of  28°  C,  the 
saturation  concentration  for  oxygen  is  8.0  mg/l.  At  a  depth  of  34  ft, 
the  hydrostatic  pressure  is  approximately  twice  that  at  the  impoundment 
surface.  Thus,  the  partial  pressure  of  oxygen  at  this  depth  is  twice 


that  at  the  surface.  Hence,  according  to  Henry's  Law,  the  saturation 
concentration  is  16.0  mg/l. 

88.  This  effect  of  hydrostatic  pressure  is  important  in  defining 
the  reaeration  process  during  venting  because  the  vented  air  travels 
with  the  discharge  water  downward  in  the  draft  tube  below  the  turbine. 

As  the  air  bubbles  travel  downward,  they  experience  increased  hydro¬ 
static  pressure  and,  as  a  result,  their  partial  pressures  increase.  At 
this  thermodynamic  state,  the  saturation  deficit,  which  is  a  measure  of 
the  force  driving  oxygen  transfer,  is  larger  than  at  atmospheric  pres¬ 
sure.  Since  the  deficit  is  larger,  more  oxygen  can  be  transferred  to 
the  water.  Thus,  increased  hydrostatic  pressure  on  air  bubbles  improves 
the  oxygen  transfer,  and  the  mathematical  description  of  the  reaeration 
process  due  to  vented  air  must  include  these  hydraulic  conditions. 

89.  Buck,  Miller,  ard  Sheppard  (1980)  developed  the  concept  of  a 
"pressure-time  history"  ot  flow  to  account  for  the  changes  in  hydro¬ 
static  pressure  as  flow  passes  through  the  draft  tube.  To  develop  such 
a  history,  a  representative  water  flow  rate  for  the  turbine  is  chosen. 
Bernoulli's  equation  (Equation  11)  and  the  continuity  equation 
(Equation  16)  are  then  applied  to  the  flow  through  the  draft  tube  to 
compute  pressures  at  several  locations  and  travel  times  between  them. 
Using  these  computations  and  assuming  atmospheric  pressure  exists  at  the 
venting  port  on  the  turbine  hub,  the  time  of  travel  for  an  average  water 
particle  can  be  plotted  against  the  pressure  which  it  experiences.  This 
is  termed  the  "pressure-time  history."  Figure  16  shows  the  pressure¬ 
time  history  for  Clarks  Hill. 

90.  While  it  is  realized  that  the  t ime-of-travel  of  a  water  par¬ 
ticle  through  a  draft  tube  is  flow  rate  dependent,  the  magnitudes  of  the 
hydrostatic  pressures  are  position  dependent  only;  thus,  the  time  value 
for  the  history  can  be  linearly  scaled  according  to  the  actual  turbine 
discharge.  For  a  discharge  that  is  smaller  than  the  one  selected  for 
history  development,  the  time  values  would  be  adjusted  with 


where 


ta  =  adjusted  time  of  flow  for  discharge  Qa  ,  sec 

t_  =  time  of  flow  for  selected  discharge  Q  ,  sec 

Qg  =  selected  discharge  for  pressure-time  history 
development,  cfs 

Q„  =  actual  turbine  discharge,  cfs 

a 

This  adjustment  to  the  pressure-time  history  allows  its  application  to 
the  range  of  turbine  operation. 


Figure  16.  Pressure-time  history  for  Clarks  Hill 


91.  A  numerical  model  was  developed  to  track  a  water  particle  as 
it  moves  through  the  draft  tube.  It  temporally  steps  through  the  pres¬ 
sure  history  of  the  draft  tube  with  very  small  time  intervals  and  solves 
Equation  19  for  final  DO  deficits  at  each  step,  accounting  for  the  im¬ 
pact  of  increased  pressure  on  the  deficit.  In  the  model,  the  final  de¬ 
ficit  for  the  previous  time  step  becomes  the  initial  deficit  for  the 
current  time  step.  In  this  manner,  the  initial  DO  deficit  is  "stepped" 
through  the  draft  tube  in  finite  increments.  At  each  of  these  time 
steps  an  improvement  of  the  release  DO,  due  to  venting,  is  achieved. 

This  stepwise  reaeration  continues  until  the  air  bubbles  reach  the  sur¬ 
face  in  the  tailrace. 

92.  The  EDM  (Equation  10)  is  then  applied  to  the  DO  deficit  for 


the  prediction  of  final  release  DO  concentration.  The  combination  of 
che  pressure-time  model  (Equation  19),  the  numerical  (computer)  tech¬ 
nique  of  stepping  through  the  pressure-time  history,  and  application  of 
the  EDM  comprises  the  numerical  model  of  turbine  reaeration  called 
VENTING.  A  simplified  schematic  of  the  computer  model  is  shown  in 
Figure  17. 

93.  To  use  VENTING,  an  estimate  of  a  gas  transfer  coefficient, 
which  is  the  unknown  quantity  in  Equation  19,  must  be  developed.  A 
least  squares  regression  analysis  of  part  of  the  available  observed  DO 
and  Q-air/Q-water  ratio  and  the  temperature  correction  (Churchhill, 
Elmore,  and  Buckingham  1962;  Tsivoglou  and  Wallace  1972)  relationship 

aT  =  a20(1.024)(T_20)  (21) 

where  a^.  and  a2Q  represent  gas  transfer  coefficients  for  water  tem¬ 
peratures  of  T  (ambient)  and  20°  C,  respectively,  were  used  to  estimate 
the  gas  transfer  coefficient.  However,  in  performing  this  analysis,  it 
was  necessary  to  include  the  turbulent  reaeration  model  (Equation  15) 
and  the  effect  of  the  pressure-time  history  on  the  deficit  ratio.  Thus, 
determining  the  estimate  of  a20  was  an  iterative  process  that  used  the 
numerical  model  VENTING. 

94.  Regression  analysis  with  the  data  taken  at  Clarks  Hill  in 
1982  resulted  in  an  estimate  of  a2Q  =  0.33/sec  .  Figure  18  shows  the 
release  DO  predicted  with  the  model  VENTING  versus  the  observed  release 
DO  concentrations  for  the  1982  field  studies.  The  standard  error  of 
estimate  for  these  data  was  0.5  mg/l.  The  maximum  prediction  error 
was  -1.1  mg/l.  Predictions  were  then  made  with  the  model  for  the  1981 
data.  These  data  were  not  used  in  the  model  development.  Figure  19 
shows  the  predicted  versus  observed  release  DO  concentrations  for  the 
1981  field  data  (also  shown  in  Tables  A8-A9).  The  standard  error  of 
estimate  for  these  predictions  was  0.3  mg/l.  The  maximum  prediction 
error  was  -0.8  mg/l.  These  results  indicate  that  the  reaeration  proc¬ 
esses  at  work  in  the  turbine  venting  can  be  described  quite  accurately 
for  the  Clarks  Hill  project. 


if 

5VJ 

*  ’J 


v 

u 

*2 

i, 

»• 


ITERATE  THROUGH 
PRESSURE  HISTORY 


INPUT  PENSTOCK  DEFICIT 


CALCULATE  NEW  DEFICIT 


ADJUST  DEFICIT  FOR 
HYDROSTATIC  PRESSURE 


CALCULATE  DEFICIT  AFTER 
TURBULENT  REAERATION 


CALCULATE  RELEASE  DO 


Figure  17.  Simplified  schematic  of  the  reaeration  computations 
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95.  The  numerical  model,  coupled  with  the  relationship  between 
efficiency  loss  and  Q-air/Q-water  shown  in  Figures  11  and  12,  allows  us 
to  estimate  the  costs  for  improving  release  DO  with  turbine  venting. 

For  this  example,  two  Q-air/Q-water  ratios  (1.5  and  3.5  percent),  four 
penstock  (upstream)  DO  concentrations,  and  two  turbine  discharges 
(3,000  and  4,500  cfs)  will  be  used  to  show  the  utility  of  the  numerical 
model  and  efficiency  loss  relationships.  The  pressure-time  history  of 
the  Clarks  Hill  project  will  also  be  used  as  input  to  the  model.  A 
complete  listing  of  the  model  and  the  input  data  to  VENTING  for  this 
example  is  given  in  Appendixes  C  and  D.  Assuming  a  water  temperature 
of  28°  C,  Table  3  displays  what  could  be  expected  at  Clarks  Hill  for  re 
lease  DO  (since  the  hydraulic  data  are  Clarks  Hill  data)  with  turbine 
venting  for  the  data  outlined  above. 


Table  3 

Predicted  Release  DO,  mg /it. 

(Water  Temperature,  28°  C;  Cs  =  8.0  mg/s.) 


Penstock  DO 
mg/ 1 

Q-Water : 

3,000  cfs 

Q-Water : 

4,500  cfs 

VV  -'•* 

VQW  *3.5* 

Qa/Qw  =1.5* 

QA/QW  =  3’5% 

0.1 

1.9 

3.4 

2.1 

3.3 

0.5 

2.2 

3.7 

2.4 

3.5 

1 .0 

2.6 

4.1 

2.8 

3.9 

4.0 

5.1 

6.2 

5.2 

6.0 

*  Q^/Qy  =  air  flow-to-water  flow  ratio. 


96.  The  improvement  in  release  DO  with  the  1.5-percent  air  flow 
is  somewhat  limited  compared  to  3.5-percent  air  flow.  The  largest  DO 
uptake  experienced  with  the  1.5-percent  air  flow  was  2.0  mg/i,  whereas 
the  maximum  DO  uptake  for  the  3.5-percent  air  flow  was  3-3  mg/5. .  Sig¬ 
nificant  improvement  for  release  DO  can  be  achieved  with  the  higher  air 
flow  rate.  However,  a  price,  in  terms  of  greater  efficiency  loss  (a 
loss  of  power  production)  and  thereby  a  loss  in  revenue,  must  be  paid. 


97.  From  Figures  11  and  12,  the  efficiency  loss  at  a  Q-air/- 
Q-water  ratio  of  1.5  percent  would  be  about  1.0  percent;  for  a  Q-air/- 
Q-water  ratio  of  3.5  percent,  the  efficiency  loss  would  be  approxi¬ 
mately  2.5  percent.  These  losses  are  a  result  of  the  air  that  was  in¬ 
troduced  into  the  flow.  Thus,  if  deflectors  were  used  to  enhance  the 
air  flow,  an  additional  loss  would  have  to  be  included  to  accurately 
assess  the  economics  of  such  a  system.  By  considering  plant  size  and 
operation,  this  efficiency  loss  can  be  converted  to  power  loss. 

98.  For  the  purpose  of  this  example,  the  effect  of  the  effi¬ 
ciency  loss  was  computed  in  terms  of  power  loss  using 

PL  =  1 . 356 ( 10~6)  eLQYH  (22) 

where 

P^  =  power  loss,  MW 

1.356(10"^)  =  constant  for  conversion  of  ft-lb/sec  to  MW 
e^  =  efficiency  loss,  percent 
H  =  gross  head,  assumed  (for  example,  146  ft) 

99.  Using  Equation  22  and  the  1.0-  and  2.5-percent  efficiency 
losses,  power  losses  due  to  venting  (losses  due  to  hub  deflectors  not 
included)  are  displayed  in  Table  4. 

Table  4 


Power  Loss  Due  to  Venting,  MW 

Efficiency  Los 

Discharge,  cfs 

% 

3,000 

4,500 

1  .0 

0.37 

0.56 

2.5 

0.93 

1.39 

100.  To  determine  the  costs  of  these  power  losses,  the  duration 
of  turbine  operation  must  be  included.  For  every  hour  that  the  turbine 
operates  at  4,500  cfs,  about  1.4  MWhr  of  energy  would  be  lost  due  to 
venting  air  into  the  turbine  at  a  rate  of  3.5  percent  of  the  water 
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flow  rate.  The  unit  price  of  megawatt-hours  could  then  be  used  to 
compute  the  actual  revenue  losses.  Additional  information  on  cost  is 
available  in  Lewis  and  Bohac  (1984). 


PART  VI:  SUMMARY  AND  CONCLUSIONS 


101.  Hydropower  is  one  of  the  cleanest  domestic  sources  of  en¬ 
ergy.  Its  use  has  been  expanding  and  will  inevitably  continue  to  ex¬ 
pand  as  nonrenewable  sources  of  energy  (fossil  fuel)  diminish  or  become 
more  expensive.  Even  though  hydropower  possesses  many  very  attractive 
attributes,  there  are  some  potential  adverse  impacts  that  can  affect 
the  quality  of  water  both  in  the  reservoir  and  downstream.  The  most 
frequently  cited  potential  adverse  impact  is  the  release  of  water  with 

a  low  DO  concentration.  This  may  occur  if  the  reservoir  is  thermally  or 
chemically  stratified. 

102.  Naturally  occurring  chemical  and  biological  processes  reduce 
the  level  of  DO  in  the  lower  levels  of  the  lake  and,  because  of  the 
thermal  stratification  (density  stratification),  this  oxygen  cannot  be 
replenished  by  reaeration  at  the  reservoir's  surface.  Hence,  the  lower 
levels  of  the  reservoir  may  become  low  in  DO  or  even  anoxic.  As  a  re¬ 
sult,  the  releases  from  a  hydropower  project  (since  hydropower  projects 
usually  withdraw  water  from  deep  in  the  upstream  pool)  can  be  very  low 
in  DO.  With  these  conditions,  the  quality  of  water  released  from  the 
hydropower  project  in  some  instances  may  be  unacceptable  relative  to 
objectives  for  the  downstream  environment  unless  measures  are  undertaken 
to  improve  the  DO  in  the  release  water. 

103.  Several  techniques  to  improve  the  quality  of  turbine  re¬ 
leases  are  available.  These  techniques  are  of  three  general  cate¬ 
gories:  forebay  systems,  tailwater  systems,  and  in-structure  systems. 
The  appropriate  reaeration  alternative  must  be  determined  on  a  case-by¬ 
case  basis.  Applicability  of  tailwater  systems  appears  limited  because 
of  the  usually  large  discharge  rate  of  hydropower  projects.  Many  of 
the  forebay  systems  partially  or  totally  break  up  the  thermal  strati¬ 
fication  in  the  reservoir.  This  may  be  unacceptable  due  to  a  possible 
change  in  the  reLease  temperature  or  mixing  of  the  distinct  water 
quality  layers  within  the  reservoir.  HypoLimnetic  aeration  or  oxy¬ 
genation  systems,  ; f  pr-  oerly  designed  and  operated,  could  provide 
improved  release  DO  . .tnout  significantly  affecting  the  reservoir's 


stratification  patterns.  However,  the  potential  for  nitrogen  super¬ 
saturation  exists  for  the  aeration  systems.  The  addition  or  retrofit 
of  selective  withdrawal  facilities  to  a  project  can  provide  the  flexi¬ 
bility  to  improve  the  release  DO;  however,  release  temperature  is 
usually  increased,  which  may  be  an  unacceptable  consequence.  Further, 
the  addition  of  a  selective  withdrawal  system  to  accommodate  the 
magnitude  of  hydropower  discharges  would,  in  most  cases,  be  very 
expensive. 

104.  Of  all  the  alternatives,  the  in-structure  techniques,  which 
usually  involve  the  injection  or  aspiration  of  air  or  oxygen  into  the 
release  flow  at  different  locations  in  the  structure,  are  the  most  at¬ 
tractive  for  release  DO  improvement.  In  particular,  for  CE  projects, 
turbine  venting  appears  the  most  applicable.  Briefly,  the  advantages 
of  a  turbine  venting  system  are: 

a.  Normally,  all  project  releases  pass  through  the  turbine, 
thereby  allowing  the  entire  project  outflow  to  be  enhanced  at 
one  location. 

b.  In  some  cases,  no  mechanical  means  or  external  power  sources 
are  required. 

c.  Turbine  venting  usually  has  no  detrimental  aesthetic  impact  on 
the  reservoir  or  tail  race. 

d.  Costs  for  a  venting  system  (capital  and  operational)  are 
usually  less  than  for  other  alternatives. 

105.  There  are,  however,  some  disadvantages  and  limitations  for 
a  turbine  venting  system,  including: 

a.  Generally,  reductions  in  turbine  efficiency  and  capacity  have 
been  observed  at  test  sites. 

b.  The  amount  of  reaeration  may  be  limited,  possibly  as  a  result 
of  the  hydraulics  of  the  release  system  of  a  hydropower 
project. 

106.  Several  alternatives  exist  with  regard  to  the  method  by 
which  a  turbine  is  vented.  For  Francis  turbines,  low  pressure  in  the 
draft  tube  at  low  operating  levels  (low  wicket  gate  settings)  causes  air 
to  be  aspirated  into  the  flow.  Ordinarily,  the  vacuum-breaker  system 
(pipe  and  valve  network  designed  to  vent  the  turbine  to  alleviate  the 
low  pressure  and  prevent  cavitation)  conducts  the  air  flow  into  the 


turbine.  However,  the  capacity  of  the  vacuum-breaker  system  is 
limited.  Significant  improvement  of  the  venting  capability  (increased 
air  flow)  can  be  achieved  by  installing  a  bell-mouth  intake  on  the  air 
supply  to  the  turbine,  thereby  avoiding  the  large  aerodynamic  losses  of 
air  flow  through  the  vacuum-breaker  system.  Even  more  air  flow  can  be 
vented  into  the  turbine  if  a  compressor  is  employed  for  forced-air  in¬ 
jection  (although  in  the  study  reported  herein,  the  compressor  was  too 
small  for  significant  improvements).  Very  large  flow  rates  of  air  may 
be  introduced  to  the  turbine  flow  if  deflector  plates  are  installed  on 
the  hub  of  the  turbine  upstream  of  the  venting  ports. 

107.  The  study  of  turbine  venting  at  the  Clarks  Hill  hydropower 
facility  was  designed  to  improve  our  understanding  of  the  reaeration 
process  during  turbine  venting  and  to  provide  guidance  on  the  appli¬ 
cation  of  this  technique.  The  results  of  the  study  indicate  that  tur¬ 
bine  venting  can  be  an  excellent  method  of  improving  concentration  of  DO 
in  releases  from  hydropower  projects.  The  following  is  a  general  list 
of  conclusions  resulting  from  the  turbine  venting  study  at  the  Clarks 
Hill  project: 

a.  Significant  reaeration  and  improvement  of  release  DO  can  be 
achieved  by  employing  turbine  venting.  In  terms  of  the  oxygen 
deficit,  this  study  indicated  that  about  30  percent  of  the 
upstream  deficit  could  be  satisfied  at  Clarks  Hill  Reservoir. 

b.  Two  processes  were  identified  with  the  uptake  of  oxygen  in 
hydropower  releases:  (1)  turbulent  reaeration  (oxygen  trans¬ 
fer  from  the  atmosphere)  in  the  tailrace  just  downstream  of 
the  draft  tube  outlet  and  (2)  aspiration  or  injection  of  air 
at  the  turbine,  resulting  in  DO  transfer  to  the  water  from  the 
air  bubbles  as  they  travel  with  the  release  through  the  draft 
tube.  From  this  study,  the  turbulent  reaeration  in  the  tail- 
race  accounted  for  up  to  6  percent  of  the  reduction  in  the  DO 
deficit . 

c.  The  effectiveness  of  the  turbine  venting  systems  studied  to 

improve  the  release  DO  was  highly  dependent  on  the  air  flow 
rate  relative  to  the  turbine  discharge  rate.  A  reduction  in 
the  upstream  DO  deficit  of  about  30  percent  was  achieved  with 
an  air  flow  rate  of  about  3  percent  of  the  turbine  discharge 

rate.  Less  reaeration  was  observed  with  lower  air  flow  rates. 

d.  The  cost  of  achieving  reductions  in  the  DO  deficit  was  re¬ 

flected  in  r  duced  turbine  operating  efficiency.  For  the 
example  gi  en  above,  an  air  flow  rate  of  3  percent  of  turbine 


discharge  resuLted  in  an  efficiency  reduction  of  approximately 
2  percent.  The  reduction  is  due  to  the  introduction  of  air  to 
the  turbine.  If  deflector  plates  are  employed  to  enhance  air 
flow,  an  additional  efficiency  loss  would  have  to  be 
included.  This  effect  on  turbine  efficiency  was  an  actual 
loss  of  generating  potential;  hence,  it  can  be  assigned  an 
actual  cost. 

e.  At  Clarks  Hill,  for  the  turbine  without  deflector  plates,  a 
slight  increase  in  turbine  efficiency  was  observed  for  very 
small  air  flow-to-water  flow  ratios  (less  than  1  percent). 

This  slight  improvement  in  efficiency  was  probably  due  to  a 
smoother  running  turbine.  However,  the  impact  on  the  release 
DO  for  these  small  air  flow-to-water  flow  ratios  was 
minimal.  This  slight  increase  in  efficiency  could  be  used  to 
offset  all  or  part  of  the  loss  in  efficiency  that  occurs  when 
venting  for  DO  enhancement.  During  periods  of  high  DO 
release,  e.g.,  late  fall,  winter,  and  spring,  venting  the 
turbine  with  small  air  volumes  could  improve  power  production 
up  to  0.75  percent  and  partially  offset  the  power  loss  during 
the  summer  months. 

f.  To  determine  if  turbine  venting  is  a  tractable  alternative,  a 
technique  to  evaluate  the  potential  improvement  in  DO  is 
required.  Based  on  the  data  collected  in  this  study,  a 
numerical  model  was  developed  as  a  predictive  tool  for 
estimating  the  effects  of  turbine  venting  on  reducing  DO 
deficits.  The  computer  model  includes  in  its  formulation  the 
reaeration  that  occurs  in  the  tailrace  due  to  turbulence  and 
the  gas  transfer  due  to  the  transport  of  air  bubbles  through 
the  draft  tube.  The  mathematical  description  of  the  former 
process  was  based  on  the  premise  that  reaeration  in  the 
tailrace  was  a  function  of  the  energy  that  was  dissipated.  It 
was  approximated  by  the  kinetic  energy  (velocity  head)  of  the 
turbine  discharge  as  the  flow  exits  the  draft  tube.  The 
latter  gas  transfer  process  was  described  by  using  a  pressure¬ 
time  history  concept  to  account  for  the  effect  of  increased 
hydrostatic  pressure  (which  changes  the  thermodynamic  state) 
on  the  air  bubbles  as  they  move  with  the  turbine  discharge 
through  the  draft  tube.  The  combination  of  these  two  process 
descriptions  resulted  in  the  numerical  model  VENTING,  which 
will  predict  the  oxygen  improvement  due  to  turbine  venting. 

108.  It  must  be  pointed  out  that  these  data  represent  the  observed 
response  of  turbines  at  the  Clarks  Hill  project.  It  is  reasonable  to 
expect  similar  responses  from  facilities  of  similar  size,  geometry,  and 
equipment.  However,  it  must  be  recognized  that  for  significantly  dif¬ 
ferent  projects,  the  applicability  of  these  predictions  may  be  limited. 
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APPENDIX  A:  TABULATED  FIELD  STUDY  DATA 
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Data  from  Unit  2,  No  Deflectors,  100%  Air 
Bell  Mouth.  1982  Study 
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Table  A7 

Data  from  Unit  Deflectors.  1001  Air 


5; 


u_  X 

cn 

i=<  £ 


t1  Qi  D  m 

ro  >  r«  3 
cu  o  o 

I  U  fo  *' 
I>  <►. 


Qj  (Tj  x-v 

+2  t-  (J 
fD  0j  — 
3  £L 

e 

01 


C  ^  £ 

3 

O 


O  —  I 

O  -  —  I 

-♦->  CD  '-v  I 

in  •  xn  \ 

!-•—«£  I 
Oi  'W  I 


3  <P  I 

«-  O  «*-  I 

-»  — *  O  I 

<r  u.  i 


ip 

I/I  D 

o  m  -*j 

t  014- 

O  X  — 


(J  DN 
^  CO  ^ 


u*i  UO  U")  U'l  tf)  U'j  Li")  Li")  ID  IT) 

«r  <r-*r<r«r'«r‘<r-r*r^r 

r*-  r*  rv  r*  r*  r-  r>  f- 


r  4  cvj  r-j  ro  ro  ro  o  <r  o 
^i>onnroncococo 

r>.  r*.  r*.  r*.  r>-  n  n  rv  N 


in  o  o  o  o  in  ID  O  O  Li") 
-Ci-O'O'O  'O'<»'O'0vO>0 


<pooo  >0-00000 

O  O  O  O  U‘j  «r  U")  o  o  o 


r^rvr^.rs.p^r^r^.  rvr-.rv 


in  in  m  in  in  o  in  o  o  o 
<i  ii)  fJ  *r  o 

ro  ro  m  n  n  f  j  n  ro  ro  ro 


O  0  0-00  CO  C4  o  O  T4  U") 

o  in  in  ro  ro  *r  ^  in  u “■  in 


oooooooooo 


rv  rv  n  co  oo  n  h 


ininr'vo-'tNLno— «rv-*r 
^ooo>oo-o^ockO 


oo^rs^ririTH^rv.  mro 
TrooNHaicO'Caio 


>o  o  «r  r\  th  r<j  rv  cl  o 
n  ro  m  n  ^  «r  t-  ld 


QJ 

tn 

L- 

fo 

IP  1 

>o 

in 

O 

o 

ro 

'Cl 

r>i 

o 

o> 

JC 

V-  1 

uo 

•o 

in 

r>. 

CN 

o- 

C  4 

o 

CN 

00 

o 

o  > 

o 

o< 

CN 

m 

00 

o 

ro 

o- 

rv 

00 

IP 

v»-  1 

cn 

04 

ro 

ro 

ro 

•«r 

ro  — i  >o  in  o  rv  co  b")  n 
rj  cn  ro  r  4  <n  cn  cn  cn  *•+  ro 


omoinoinoinoin 

■^*<riniD'O'Or'vr>s0000 


,4's.l 


li 

ill; 


1 


fefcifc 

$& 


$ 


b-^vv 


y 

I,V<'V 

«V  V 


r;W  •  -  ‘1#  "V  V 


$ 


$ 


$ 

$ 


v 


ill 

V/i 


jt'l 


ti'J 


i**1 


V.J 

JS 


,<X*1 

•y,1 

!!v, 

V»j: 

v«»i 

■‘O' 


<  v 

4 


& 
;k 


«►»  i>  vi 
Qi  m  X 
•  IT 
O  01  £ 


CO 

L* 

OO  O 

<X  4J 

o 

0)  <D 

<— <  r— I 

X>  C*-H 

CTJ  0> 

E-  Q 


*o 

a>  ~ 

■♦>  *  — * 

U  O 
^  -  tn 

■  X7  u  fi 

0)  0i 

a<  t- 

«-  a. 


3  x? 

o  b 

Cl  >  •  *— » 

c-  O 
04  •  VI 

in  a  « 
-Q 
O 


U  —  I 

O  •  —♦  I 

O 

Ip  -  VI  l 

cam  i 

a 4  i 


«-  3  IP  I 

^  O  «•-  I 

<x  -H  O  1 


ai  (L  3  I 
3  4->  ac  i 

or»^  i 


ro  ip 

JO  <+- 

u  o 

IP  ■— 


Ip 

Ip  *o 
O  ro  +> 
<-  04  <W 
O  X  ^ 


0l  0/  I 

-X  ~  I 

U  ffl  N* 

H0  v  I 


''O-O'-O'O'-oo-O'Orv  'Orx'Orv.sorxrxrv.rvrxrx 

p  in  ui  m  in  ut  lt  Li  ip  in  in  ut  LT  m  in  in  n  in  in  ip 
rxrxrxrxrxrvr^rxrxrxrxrxrxrxrxrxrxrxr'' 


ooooooooo  ooooooooooo 
OO'OOOOO'— *rooroo  ^  o  o«  O'  O'  O'  O'O' 

m  oo  O'  O'  O'  O'  00  O'  nn  o*  CO  O'  QQ  O'  00  fp  nr>  m  pn  tp 


u^rx'OOfori'O'OO'  ororocoo'o-ooo  n  rg 
«rro^orx-or>i— <o>  — *o— <o—irN4— «-HOO'ao 


oooinmomoao  oomoinoooooo 

iS'OroH00'O'Obim'O^N(M)'<r  noO'.O  n 


oooooooooooooooooooo 

O'  O'  O'  O'  O'  O'  O'  O'  O'  O'  O'  O'  O'  O'  O'  O'  O'  O'  O'  O' 

<n  r*  cm  cm  cm  r> i  cm  cm  cm  cm  cm  cm  cm  cm  cm  cm  cm  cm  cm  cm 


n  h  n  o  n  n  h  r  ^  rjoioo^-r-xor^oinro 

O'— lODO'^rcDaoroO'  O'— «o^'oo^rxO'C^oo 
iP'0<rMriHininn^«fp^inT  ▼  n  n  ^ 


'OO'O'^- i«OHO'M^ooDHooNcomnin  o  o 
r  <  m  o  rj  v  h  h  n  'O  in  m 

—iin-v^r^’^rrvo'— iroT'O'orxr^rxrxrxfx.rx 
nroMfowromrot  *r  v  <r  ^  ^  v  v 


n  h  o  h  ih  n  'O  — «  -<r  co  o  n  rv.'OO  o  ro  o  w  n 
©  h  n  n  rj  m  w  y-j  v  r^ro^oo^ao  -©O'O'  O'  o* 

Hin^V^^N^HrtliTNCDOO^HHHH^ 

nnroMnwnnT  ^  -*r  ^pppinininin 


HP^OrnM^OOI/)N©(MO(MP^«(NHp'M 

^-(^nON^OnOWOnOCOOOOODNN 

*0'0'0«o>oin'0'0<‘  <r  ^ 

mwrowwfOMnn  nw  w  wrow  fomro  ro  ro 


oinmininmotDomoinoinoooooo 

goininminiD'O'or^rvQoaiO'O'OOOOOo 


& 


w,» 


«  j,"  * 

g& 


Table  A9 

Data  from  Unit  4.  Deflectors.  01  Air 


fSi 

'ill 

i 


yi 


•/»*< 

$ 


+>  3  xn 
a i  oi  i 
m  ai 
O  Oj  C 


r» 

t- 

l  (0  ^ 

<-  u 
m  (u 
3  a. 


X3 

Oj  — 

u  o  s 

•  tn 

«  ~o  i=.  m 

ft»  a»  — 

at  u 


*  10 

O  at  -N 
C  >  •  -• 
«-  O  ''v 
«l  •  w 
(A  O  fl 
JO  — 

o 


u  » 

O  -  — •  l 

O  \  I 

at  •  tn  i 

C  O  I  f 

Of  —  I 

a.  i 


t  J  IA  > 

M  O  %-  I 

C  ^  O  I 


u  o  ~+ 

tl  0.  3 
J  « 

0  3  — 
CL  O 


(0  at 

-c  %- 

u  u 

at  — 


at  *o  ~ 
o  m  ♦> 
a»  v 
o  x  — 


•i  * 

-L 

OWN 
CD  — 


-o  'O'O  'O,o>or,s'orv  rs  rs.  *o  'O 

in  inmmaiifjmminmuiinin 
^rsrvrvrvrsNrvr\rsNrNrs 


ooooooooooooo 
Os  O  O'  HODChrOOMO^On 

CD  OK  GO  O'  00  00  00  O'  00  O'  00  0>  CD 


moooomooooooo 


m  po  m  to  ro 


m  m  m  fo  ro  ro 


ooooooooooooo 

O'  OO'O'O'OO'O'O'OO'O'O' 

ra  rj  cm  cn  p>i  ra  ra  ra  ra  ra  ra  cm  ra 


^intNinw^oo^oD^Naoro 

f^NN^TfO^vfOO^OOin 

(Nfinaocoo^rof-«'ONN 


O'O^rv^rv.uTO'm'Oinom 

'Ocoh  onm  o 
ra«rmrvrvO'-*ro^rsO*o-« 
rororororofOT^'T  t  t  ^  in 


rv>o*^o^ps>oropNrapN^4(>s 

oo«-*ooorooroo^om 

rtMMWIOMMmoniOMrt 


ommooinoinoiDomo 

uiinm'O'O'ONNODoooo 


,V(Y 

M 

Fife 

ft;' 

Mh, 


m 

moV 

,w» 

B 


0 


m 


v 
i:n\ 


pv'V 

$S 

m 


|§: 


i 

♦V  <  • 


•<.  S 


APPENDIX  B:  GRAPHICAL  PRESENTATION  AND  ANALYSIS 
OF  FIELD  DATA 
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Wicket  Gate  Setting,  7. 


Figure  B1.  Power  output  versus  wicket  gate,  power  output 
adjusted  to  146-ft  gross  head,  Unit  2,  no  air  flow 


Wicket  Gate  Setting,  7 


Figure  B2.  Discharge  versus  wicket  gate,  discharge 
adjusted  to  146-ft  gross  head,  Unit  2,  no  air  flow 
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Wicket  Gete  Setting,  X 


igure  B3.  Power  output  versus  wicket  gate,  power  output  adjusted 
to  146-ft  gross  head,  Unit  2,  vacuum  breaker  open 
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Figure  B4.  Discharge  versus  wicket  gate,  Discharge  adjusted 
to  146-ft  gross  head,  Unit  2,  vacuum  breaker  open 
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Figure  B5.  Power  output  versus  wicket  gate,  power  output  adjusted 
to  146-ft  gross  head,  Unit  2,  bell-mouth  intake  air  supply 
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Figure  B6.  Discharge  versus  wicket  gate,  discharge  adjusted 
to  146-ft  gross  head,  Unit  2,  bell-mouth  intake  air  supply 


Figure  B7.  Power  output  versus  wicket  gate,  power  output  adjusted 
to  146-ft  gross  head,  Unit  2,  blower-supplied  air 
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Figure  B 1 3 .  Power  versus  wicket  gate,  power  output  adjusted 
to  146-ft  gross  head,  Unit  4,  bell-mouth  intake  air  supply 


Figure  B14 .  Discharge  versus  wicket  gate,  discharge  adusted 
to  146-ft  gross  head,  Unit  4,  bell-mouth  air  supply 
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Figure  B15.  Efficiency  versus  wicket  gate,  Unit  2 
.  =  no  air  flow,  o  =  vacuum  breaker  open 
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Figure  B 1 6 .  Efficiency  versus  wicket  gate,  Unit  2 
.  -  no  air  flow,  o  =  bell-mouth  intake  air  supply 
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Figure  B 1 9 -  Efficiency  versus  wicket  gate,  Unit  4 
.  =  no  air  flow,  o  =  bell-mouth  intake  air  supply 
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Figure  B20.  Discharge  versus  efficiency,  Unit  2 
.  =  no  air  fLow,  o  =  vacuum  breaker  open 
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Figure  B25.  Water  flow  and  air  flow  versus  wicket  gate  setting, 
Unit  2,  vacuum  breaker  open,  *  =  air  flow,  o  =  water  flow 
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Figure  B26.  Water  flow  and  air  flow  versus  wicket 
gate  setting,  Unit  2,  bell-mouth  intake  air  supply 
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Figure  B27.  Water  flow  and  air  flow  versus  wicket  gate  setting, 
Unit  2,  blower-supplied  air,  *  =  air  flow,  o  =  water  flow 
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Figure  B28.  Water  flow  and  air  flow  versus  wicket  gate  setting, 
Unit  4,  vacuum  breaker  open,  *  =  air  flow,  o  =  water  flow 
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Figure  B29.  Water  flow  and  air  flow  versus  wicket 
gate  setting,  Unit  4,  bell-mouth  intake  air  supply, 
*  =  air  flow,  o  =  water  flow 
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PROGRAM  TURAER< INPUT. OUTPUT, TAPES. TAPE6=0UTPUT) 


1000 
1001* 

1002*  APRIL  1983  WES-HS3 

1003*  purpose: 

1004*  PREDICTION  OF  THE  DISSOLVED  OXYGEN  UPTAKE  OF 

1005*  DRAFT  TUBE  HYDROELECRTRIC  TURBINE  AERATION  SYSTEMS. 

1006******************************************************************* 

1007* 

1008* 

1009********INPUT. 

1010* 

1011*  QAIR  -  AIR  FLOW  RATE  (FT3/SEC ) 

1012*  QWATER  -  WATER  FLOW  (TURBINE  DISCHARGE)  RATE  (FT3/SEC ) 

1013*  TEMP  -  OUTFLOW  WATER  TEMPERATURE  (DEG  C) 

1014*  DOICON  -  INITIAL  (UPTAKE)  DISSOLVED  OXYGEN  CONCENTRATION  (PPM) 

1015*  ELCL  -  CENTERLINE  ELEVATION  OF  DRAFT  TUBE  OUTLET 

1016*  BETA  -  ENERGY  DISSIPATION  COEFFICIENT  FOR  TURBULENCE 

1017*  ALPHA  -  GAS  TRANSFER  RATE  COEFFICIENT  FOR  VENTING 

1018*  AREA  -  AREA  OF  DRAFT  TUBE  OUTLET 

1019*  NTEST  -  TEST  NUMBER 

1020*  DBASE  -  FLOWRATE  FOR  WHICH  TIME-PRESSURE  HISTORY 

1021*  WAS  DEVELOPED. 

1022*  NPOINTS  -  NUMBER  OF  TIME-HISTORY  POINTS  (PAIRS) 

1023*  ELTW  -  TAILWATER  ELEVATION  (FT) 

1024*  •  FLOTIME  -  AIR-WATER  CONTACT  TIME  (SEC) 

1025*  BELT  -  TIME  STEP  USED  IN  SUMMATION  OVER  THE  TOTAL 

1026*  AIR-WATER  CONTACT  TIME  (SEC) 

1027*  PRESSTM  -  TIME  HISTORY. IE.  DRAFT  TUBE  PRESSURE(PSI )  VERSES 

1028*  TIME(SEC) 

1029*  DOSCON  -  SATURATED  DISSOLVED  OXYGEN  CONCENTRATION  (PPM) 

1030* 

1031*************************************************************** 

1032* 

1033*******0UTPUT. 

1034* 

1035*  DOFCON  -  FINAL  (RELEASE)  DISSOLVED  OXYGEN-  PREDICTED- 

1036*  CONCENTRATION  (PPM) 

1037* 

1038******************************************************* *********** 

1039* 

1040*********PARAMETERS. 

1041* 

1042  PARAMETER  NSTEP  =  1000.  NPRESS  -  10) 

1043* 

1044*  NSTEP  -  NUMBER  OF  PRESSURE  POINTS  OVER  THE 

1045*  TOTAL  AIR-WATEP  CONTACT  TIME. 

1046*  NPRESS  -  NUMBER  OF  PAIRS  OF  INPUT  PRESSURE-TIME 

1047*  POINTS. 

1048************** ******  ***** «******«« *(«*«<***** ******** ********** 

1049* 

1050  COMMON  /CONTROL/  PRESS ( NSTEP  '  •  PRESSTM (NPRESS . 2  • 


1052 

COMMON  /  AA  /  IFILE >  OFILE.  END.  NPOINTSf 

NTEST .  QBASE 

1053 

COMMON  /  B6  /  ALPHA.  RATIO.  QAIR.  QUATER. 

DELT.  FLOTIME.  PATM 

1054 

COMMON  /  CC  /  ELCL.  ELTU.  DOICON.  DOFCON. 

TEMP 

1055 

COMMON  /  DD  /  OFIRST.  QDEBUG.  QERROR*  QSTQP.  QECHO.  QCHECK 

1056 

COMMON  /  EE  /  AREA.  BETA 

1057 

INTEGER  OFILE.END 

1058 

LOGICAL  QDEBUG.  QFIRST .  QSTOP.  QERROR.  QECHO.  QCHECK 

1059* 

1060****** 

SET  UP  PROGRAM  CONTROL  VARIABLES  ****** 

1061* 

1062 

IFILE  =  5 

1063 

OFILE  =  6 

1064 

QFIRST  =  .TRUE. 

1065 

QDEBUG  =  .FALSE. 

1066 

QSTOP  =  .FALSE. 

1067 

QERROR  = 

.FALSE. 

1068 

QECHO  = 

.FALSE. 

1069 

QCHECK  = 

.TRUE. 

1070* 

i07i************************************************* 
1072*******  READ  INPUT  FROM  DATA  FILE 


1073* 

1074  50 

CALL  TREAD 

1075* 

1076 

IF  (QERROR)  STOP 

1077* 

1078 

CALL  HISTORY 

1079* 

1080 

CALL  AERATE 

1081* 

1082 

CALL  FINALC 

1083* 

1084***CYCLE 

BACK  TO  READ  ABDITIANAL  DATA  SETS 

1085* 

1086 

WRITE  (  OFILE.  200  ) 

1087 

GO  TO  50 

1088* 

1089 

1090 


200 


FORMAT  <  1H 
END 


1092 

1093* 

1094* 

1095* 

1096 

1097 
1098* 

1099 

1100 
1101 
1102 

1103 

1104 

1105 

1106 
1107 
1108* 

1109 

1110 
1111 
1112 

1113 

1114 

1115 

1116 

1117 

1118 

1119 

1120 


SUBROUTINE  TREAD 


READ  THE  INPUT  DATA  FOR  TURBINE  AERATION  SYSTEM 


PARAMETER  (  NSTEP  =  1000. 
IMPLICIT  INTEGER  (X) 


NPRESS  = 


COMMON  /  CONTROL  i  PRESS ( NSTEP ) .  PRESSTM ( NPRESS . 
COMMON  /  AA  /  IFILE.  OFILE.  END.  NPOINTS.  NTEST. 
COMMON  /  BB  /  ALPHA.  RATIO.  GAIR .  QWATER.  BELT, 
COMMON  /  CC  /  ELCL.  ELTU.  DOICON.  DOFCON.  TEMP 
COMMON  /  DD  /  QFIRST ,  QDEBUG.  QERROR »  QSTOP.  GEC 
COMMON  /  EE  /  AREA.  BETA 
DIMENSION  TITLE  (7),  DUMMY  (20) 

INTEGER  OF ILE. END. DUMMY. CHECK. CHECK1. TITLE 
LOGICAL  QFIRST. QDEBUG. QSTOP. QERROR. QECHO.  QCHECK 


IFILE.  OFILE.  END. 
ALPHA.  RATIO.  GAIR. 
ELCL.  ELTU.  DOICON. 


PRESSTM(NPRESS»2) 

NPOINTS.  NTEST.  QBASE 
QWATER.  BELT,  FLOTIME.  PATH 
DOFCON.  TEMP 


/  QFIRST.  QDEBUG.  QERROR.  QSTOP.  QECHU.  QCHECK 


XENGL. 
XSTOP. 
XDATA. 
XMETR » 
XDELT » 
XBLANK » 
XTIME. 
XWATER . 
XQAIR  , 
XTAIL . 
XBETA, 
XDOI  / 


XF ILE  /'4HENGL.  4HFILE  / 
XDEBUG  /4HST0P ,  4HDEBU  / 
XPRIN  /  4HDATA.  4HPRIN  / 
XUNITS  /  4HMETR .  4HUNIT  ./ 
XINTER  /  4HDELT.  4HINTE  / 
XALPHA  /  4HBLAN.  4HALPH  / 
XTEMP  /  4HTIME.  4HTEMP  / 
XBASE  /  4HUATE.  4HBASE  / 
XF'RESS  /  4HAIR  .  4HPRES  / 
XELCL  /  4HTAIL.  4HCENT  / 

X AREA  /  4HBETA.  4HAREA  / 
4HD0  I  / 


1121  ********************************************************** 


1122* 

1123* 

1124 

1125 
1126* 
1127* 
1123* 

1129 

1130 

1131 

1132 


ADDITIONAL  DATA  SETS  IF  QFIRST 


.FALSE. 


(QSTOP.  STOP 
•:  .NOT.  QFIRST 


GO  TO 


READ  TITLE  AND  FILE 


CONTROL  DATA 


IF  (CHECK 
READ  •  IF  I l 


aFILE-  505.  TITLE 
( IF  II  E .  510)  CHECK.  I  -ILE. 
'.IFILE.  505)  CHECK 
HECK  ,EQ.  XDEBUG)  QDEBUG  = 
•  IFILE.  505)  CHECK  t 


OFILE 


.TRUE. 


1134 

1135 
1136* 
1137* 
1133* 

1139 

1140 

1141 

1142 

1143 


CHECK!  . NE . 
CALL  ERROR 


'ENGL  .AND.  CHECK1  . NE  , 
CHECK.  XUNITS' 


nMETF: . 


ECHO  RP!NT 


READ 
QECHC 
IF  ■  . 

QECHC 
REWIND 


IFIlE.  505 
-  CHECK  .EG. 
NCT.  QECHO  . 
=  .FAlSE. 
IFILE 


CHECK 
XPRIN 
GO  TO  140 


I 


1144 

LINE  =  1000 

1145 

WRITE  (  OFILE,  600  ) 

1146  100 

CONTINUE 

1147 

READ ' (  IFILE,  505,  END  -  110  )  Dt 

1148 

WRITE  (  OFILE,  610  )  LINE,  DUMMY 

1149 

LINE  =  LINE  +10 

1150 

GO  TO  100 

1151  110 

CONTINUE 

1152 

REWIND  IFILE 

1153 

DO  130  I  =  1,  5 

1154 

READ  (  IFILE,  505  )  DUMMY 

1155  130 

CONTINUE 

1156  140 

CONTINUE 

1157* 

1158*  SET 

ATMOSPHERIC  PRESSURE  (UNITS) 

1159* 

1160 

PATH  =  1.03323 

1161  IF  (CHECKl  .EG.  XENGL)  PATH  =  14.695? 


1162* 

1163  READ  (IFILE,  525)  CHECK.  ELCL 

1164  IF  (CHECK  .NE.  XELCL )  CALL  ERROR  (CHECK.  XELCL ) 

1165  READ  (  IFILE.  525  )  CHECK.  AREA 

1166  IF  (  CHECK  .  NE .  XAREA  )  CALL  ERROR  (  CHECK.  XAREA  ) 

1167  READ  (IFILE.  520)  CHECK.  DELT 

1168  IF  (CHECK  . NE .  XINTER)  CALL  ERROR  (CHECK.  XINTER) 

1169  READ  (IFILE.  525)  CHECK.  FLOTIHE 

1170  IF  (CHECK  .NE.  XBASE)  CALL  ERROR  (CHECK,  XBASE) 

1171  READ  (IFILE,  525)  CHECK,  ALPHA 

1172  IF  (CHECK  .NE.  XALPHA )  CALL  ERROR  (CHECK,  X ALPHA) 

1173  READ  (  IFILE,  525  )  CHECK,  BETA 

1174  IF  (  CHECK  .NE.  XBETA  )  CALL  ERROR  (  CHECK,  XBETA  ) 

1175  READ  (IFILE,  505)  CHECK 
1176* 

1177*  PARAMETER  CHECK 
1178* 

117?  TSTEP  =  FLOTIME/DELT 

1180  IF  (TSTEP  .LT.  NSTEPi  GO  TO  105 

1181  WRITE  (  OFILE,  530  >  TSTEP 

1182  STOP 


1183  105  CONTINUE 

1184* 

1185*  READ  THE  NUMBER  OF  DATA  POINTS 
1186*  IN  PRESSURE-TIME  HISTORY 
1187* 

1188  IF  '.CHECK  .NE.  XF'RESS )  lALl  ERROR  (CHECK,  (PRESS 1 

1189  READ  (IFILE,  515)  CHECK,  NPOINTS 

1190  IF  .CHECK  .NE.  XTIME,  CALL  ERROR  ( CHECK • ‘TIME  • 

1191  NPOINTS  -  NPOINTS  +  1 

1192* 

1193*  PARAMETER  CHECK 

1194* 


1195 

1196 


IF  'NFO.’-'TS 
WRITE  '  Of  . 


F.  NPRESS)  GO  TO  33 
540  >  NPOINTS 


C6 


•  .  3 


1197  STOP 

1198  33  CONTINUE 

1199* 

1200*  READ  THE  PRESSURE-TINE  HISTORY 
1201* 

1202  DO  20  1=1 .NPQINTS  -  1 

1203  READ  (IFILE, 550)  PRESSTM< I , 1 ) ,PRESSTN( I , 2 ) 

1204  20  CONTINUE 

1205* 

1206  READ  <  IFILE,  525  )  CHECK»  QBASE 

1207  IF  <  CHECK  ,NE .  XBASE  )  CALL  ERROR  <  CHECK ,  XBASE  ) 
1208* 

1209*  BEGIN  DATA  THAT  CAN  BE  CHANGED 
1210*  TO  REPRESENT  ADDITIONAL  DATA  SETS. 

1211*  ALL  OF  THE  INPUT  DATA  ABOVE  IS 
1212*  UNCHANGED  FROM  ONE  DATA 
1213*  SET  TO  ANOTHER. 

1214* 

1215  1000  CONTINUE 

1216  READ  (IFILE.  560)  TITLE.  NTEST 

1217  READ  (IFILE.  525)  CHECK.  ELTW 

1218  IF  (CHECK  ,NE.  XTAIL )  CALL  ERROR  (CHECK.  XTAIL) 

1219  READ  (IFILE.  525)  CHECK,  TENP 

1220  IF  (CHECK  .NE.  XTENP )  CALL  ERROR  (CHECK,  XTEMP ) 

1221  READ  (IFILE,  520)  CHECK,  QAIR 

1222  IF  (CHECK  . NE.  XQAIR)  CALL  ERROR  (CHECK,  XQAIR) 

1223  READ  (  IFILE,  520  )  CHECK,  QHATER 

1224  IF  (  CHECK  .NE.  XWATER  )  CALL  ERROR  (  CHECK.  XWATER  ) 

1225  READ  (IFILE,  520)  CHECK,  DOICON 

1226  IF  (CHECK  .NE .  XDOI)  CALL  ERROR  (CHECK,  XDOI) 

1227  READ  (IFILE,  520)  CHECK 

1228  IF  (CHECK  .NE.  XDATA)  GO  TO  60 

1229  BACKSPACE  IFILE 

1230  QFIRST  =  .FALSE. 

1231  RETURN 

1232  60  CONTINUE 

1233  IF  (CHECK  .EQ.  XSTOP)  QSTOP  =  .TRUE. 

1234  IF  (CHECK  .NE.  XSTOP)  CALL  ERROR  (CHECK,  XSTOP) 

1235  RETURN 

1236  505  FORNAT  (  20A4  ) 

1237  510  FORMAT  (  A4,  6X,  215  ) 

1238  515  FORMAT  (  A4,  16X,  15  ) 

1239  520  FORMAT  (  A4,  6X,  F10.2  ) 

1240  525  FORMAT  (  A4,  21X,  F10.2  ) 

1241  530  FORMAT  \  'PARAMETER  NSTEP  IS  TOO  SMALL',/, 

1242  *  'SHOULD  BE  ',14.  '  RECOMPLILE  ) 

1243  540  FORMAT  ■'  ,  'PARAMETER  NPRESS  IS  TOO  SMALL  '•/, 

1244  *  'SHOULD  BE  ',  15  , 

1245  550  FORMAT ( 2F 1 0 . 2 ) 

1246  560  FORMAT  (  7A4,  A7  ) 

1247  600  FORMAT  (  1H1  ) 

1248  610  FORMAT  •  10X,  16.  7X,  3H**«,  20A4  • 

1249* 


1252  SUBROUTINE  HISTORY 

1253* 

125**  COMPUTES  THE  PRESSURE  HISTORY  OVER  THE  TOTAL  AIR-WATER 
1255*  CONTACT  TIME  AT  EACH  TIME  STEP(DELT). 

1256* 

1257  PARAMETER  (NSTEP  =  1000.  NPRESS  =  10) 

1258* 

1259  COMMON  /CONTROL/  PRESS ( NSTEP ) .  PRESSTM < NPRESS . 2 ) 

1260* 

1261  COMMON  /  AA  /  IFILE.  OFILE.  END.  NPOINTS.  NTEST.  DBASE 

1262  COMMON  /  BB  /  ALPHA,  RATIO,  QAIR,  QUATER.  BELT,  FLOTIME,  PATH 

1263  COMMON  /  CC  /  ELCL,  ELTW,  DOICON,  DOFCON,  TEMP 

1264  COMMON  /  DD  /  GFIRST,  QDEBUG,  QERROR ,  QSTOF ,  QECHO 

1265  COMMON  /  EE  /  AREA,  BETA 

1266* 

1267  LOGICAL  QDEBUG 

1268* 

1269  DIMENSION  ARRAY  ■  NPRESS  > 

1270  INTEGER  START , END , OF ILE 
1271* 

1272*  LOOP  OVER  EACH  INPUT  PRESSURE-TIME  POINT  IN  'ARRAY1 
1273*  AND  GENERATE  PRESSURE  POINTS  FOR  EACH  TIME  STEP ( DELT ) , 

1274* 

1275* 

1276*  SCALE  THE  INPUT  TIME  POINTS  BASED 
1277*  UPON  THE  RATIO  QBASE/QWATER 
1278* 

1279  NF'Ml  =  NPOINTS  -  1 

1280* 

1281  SCALE  =  QBASE  /QUA TER 

1282* 

1283  DO  10  I  =  1,  NPM 1 

1284  ARRAY  (  I  ■  =  PRESSTM  (1,1)  *  SCALE 

1285  10  CONTINUE 
1286* 

1287*  SET  LAST  TIME-HISTORY  POINT  TO  ATMOSPHERIC  USING 
1288*  RISE  VELOCITY  OF  BUBBLES  TO  BE  2.0  FT/SEC. 

1289* 

1290  VRISE  =  2.0 

1291  ARRAY  <  NPOINTS  •  =  ARRAY  (NF’Ml  )  +  (  ELTW-ElCL  > 'VRISE 

12^2  PRESSTM  (  NPOINTS,  2  )  -  PATM 

1293* 

1294  IF  ■.  ODE  BUG  -  UF  I  T  E  ■  GFILE.500  )  (  ARRA  •  I  • ,  I  -  i ,  NFO I N  TS  ■ 

1295  IF  ■  QDEBUG  <  WRITE  OFILE -505  )  !  PRESSTM'  I<  2-, 1  =  1, N  POINTS 

1296* 

1297  LOOP  =  1 

1298  PRESS'  1  )  -  P  R  E  S  S  T  M  ( 1  -  2  ■ 

1299* 

1  300  DO  9V  K  =  1 ,  N P M 1 

1301* 

1302*  FIND  THE  NUMBER  OF  TIME  STEPS  BETWEEN  TWO  INPUT 
i  303*  TIME-HISTORY  -Ol-iiS 


3 

1305 
1306* 
1307*  F 
1308* 

1309 

1310 
1311* 
1312*  F 
1313* 

1314 

1315 
1316* 
1317 
1318* 

1319 

1320  60C 

1321 

1  322** 
1323* 

1324  55 

1325 

1326  30 
1327* 

1328 

1329 

1330 
1331* 

1332  99 

1333  30C 

1334  503 

1335  5 1 ( 

1336  52C 
1337* 
1338 


JUMP  =  NINT  (  (  ARRAY ( K+ 1 )  -  ARRAY (K))  /  CELT  ) 

THE  PRESSURE  INCREMENT 

BUMP  =  (  PRESS TM(K+1 » 2 )  -  PRESSTM ( K , 2 ) )  /  FLOAT(JUMP) 
IF  (  QDEBUG  )  WRITE  (  OFILE. 510  )  JUMP.  BUMP 

IN  THE  COMPLETE  PRESSURE  ARRAY 

START  =  LOOP  +1 
END  =  LOOP  +  JUMP 

IF  (END  .LE.  NSTEP)  GO  TO  55 

WRITE  OFILE.  600  > 

FORMAT  (  .  NSTEP  IS  TOO  SMALL  -  RECOMPILE 

STOP- 


DO  80  J  =  ST ART . END 

PRESS ( J )  =  PRESS ( J-l ;  +  BUMP 

CONTINUE 

LOOP  =  LOOP  +  JUMP 

IF  <  QDEBUG  -  WRITE  (  OFILE. 520  »  START.  END. 

(  PRESS ( J 1  .  J  -  START*  END  ■ 

CONTINUE 

FORMAT  (  23H*»  SCALED  TIME  ARRAY  *».//.  10F8.2  1 
FORMAT  (  27H»*  SCALED  PRESSURE  ARRAY  10F6.2  ) 

FORMAT  <  8H JUMP  =  .  1 4 . 3X . 8HBUMP  =  .F6.2  ) 

FORMAT  (  21HFRESSURE  HISTORY  FROM . 1 4 . 2HT Q . 1 4 .  . (  10F8.2)  ) 

RETURN 


<v< 
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1341  SUBROUTINE  AERATE 

1342* 

1343*  COMPUTES  THE  DISSOLVED  OXYGEN  UPTAKE  OF 
1344*  TURBINE  REAERATION  SYSTEMS. 

1345* 

1346*  JUNE  1982  -  WES-HS3 

1347* 

1348  PARAMETER  <NSTEP  =  1000.  NPRESS  =  10) 

1349  COMMON  /CONTROL/  PRESS(NSTEP) »  PRESSTM(NPRESS»2) 

1350* 

1351  COMMON  /  AA  /  IFILE.  OFILE.  END.  NPOINTS.  NTEST.  DBASE 

1352  COMMON  /  BB  /  ALPHA.  RATIO.  QA1R.  QHATER.  BELT.  FLOTIME.  F'ATM 

1353  COMMON  /  CC  /  ELCL.  ELTW.  DOICON.  DOFCON.  TEMP 

1354  COMMON  /  DD  /  QFIRST .  QDEBUG,  QERROR.  QSTOP.  QECHO 

1355  COMMON  /  EE  /  AREA.  BETA 

1356* 

1357  INTEGER  END. OFILE 

1358  LOGICAL  QDEBUG 
1359* 

1360  RATIO  =  QAIR  /  QWATER 

1361* 

1362*  CALCULATE  THE  SATURATION  CONCENTRATION 
1363* 

1364  DOSCON  =  1.0  /  (  0.00209  *  TEMP  +  0.06719  ) 

1365* 

1366*  AIR  ASPIRATION/INJECTION  MODEL  PER  TIME  STEP 
1367* 

1368  CONSTO  =  <  ALPHA  )  *  (1  .024  **  (TEMP-20., <)  *  RATIO  *  BELT 

1369* 

1370*  CALCULATE  PENSTOCK  DEFICIT 
1371* 

1372  CS02  =  DOSCON  *  PRESS(l)  /  PATM 

1373  DI  =  CS02  -  DOICON 
1374* 

1375*  ITERATE  THROUGH  PRESSURE  HISTORY 
1376* 

1377  IF  (QDEBUG)  WRITE < OFILE . 200 ) 

1378  K  =  END  -  1 

1379  DO  100  I  =  1,  K 
1380* 

1381*  CALCULATE  NEW  DEFICIT 
1382* 

1383  DF  =  DI  *  EXP (-CONSTO) 

1384* 

1385*  ADJUST  DEFICIT  FOR  HYDROSTATIC  PRESSURE 
1386* 

1387  DOFCON  =  DOSCON  *  PRESS'D  /  PATM  -  DF 

1388* 

1389  DI  =  DOSCON  *  PRESSU  +  D  /  PATM  -  DOFCON 

1390* 

1391 

1392 


BB  =  DOSCON  *PRESS;i+l)  /  PATM 

IF  (  QDEBUG  )  WRITE  (OFILE. 50)  DI.  BB.  DOFCON.  DF 


1393* 

1394  100  CONTINUE 
1395* 

1396*  CALCULATE  VELOCITY  HEAD  AT  DRAFT  TUBE  OUTLET 
1397* 

1398  DELENG  =  <  QWATER  /  AREA  )  **  2.0  /  64.348 

1399* 

1400*  CALCULATE  DEFICIT  AFTER  TURBULENT  REAERATION 
1401* 

1402  DF  =  DI  *  EXP  (  -BETA* ( 1 .024**( TEHP  -20 . )  XDELENG) 

1403* 

1404*  CALCULATE  RELEASE  DO 
1405* 

1406  DOFCON  =  DOSCON  -  DF 


1406 

1407* 

1408 

1409* 

1410  200 

1411  50 

1412 
1413* 


RETURN 


FORMAT ( '  4X» 'DI ' >7X  > ' BOS ' * 7X» 'DOF ' »5X» ' DF ' ) 

F0RMAT(4(F8.3»2X) ) 

END 


1414  SUBROUTINE  FINALC 

1415* 

1416*  OUTPUT  FOR  -SUBROUTINE  AERATE 
1417* 

1418* 

1419  COMMON  /  AA  /  IFILE.  OFILE.  END.  NPOINTS.  NTEST »  QBASE 

1420  COMMON  /  BB  /  ALPHA.  RATIO.  QAIR.  QWATER.  BELT.  FLOTIME.  PATM 

1421  COMMON  /  CC  /  ELCL.  ELTU.  BOICON.  DOFCON.  TEMP 

1422  COMMON  /  DD  /  QFIRST .  QBEBUG.  OERROR.  QSTOP.  QECHO.  QCHECK 

1423* 

1424  INTEGER  OFILE.END 

1425  LOGICAL  QFIRST.  QCHECK 

1426  IF  (  QCHECK  )  WRITE  <  OFILE.  600  ) 

1427  IF  (  QCHECK  )  WRITE  <  OFILE.  610  ) 

1428  IF  (  QCHECK  )  WRITE  <  OFILE.  620  ) 

1429  IF  (  QCHECK  )  WRITE  (  OFILE.  622  ) 

1430  IF  (  QCHECK  )  WRITE  <  OFILE.  625  ) 

1431  QCHECK  =  .FALSE. 

1432  WRITE  <  OFILE.  630  )  NTEST. ELTW. TEMP. QWATER. QAIR. BOICON. DOFCON 

1433  600  FORMAT  (  1H1  ) 

1434  610  FORMAT  <  /// » 2X » 4HTEST , 4X » 4HTAIL » 5X . 1 1HTEMPERATURE » 

1435  *  4X.5HWATER.4X.3HAIR.4X.7HINITIAL.6X.6HFINAL  ) 

1436  620  FORMAT  <  1X.6HNUMBER.3X.5HWATER. 19X.4HFL0W.5X.4HFL0W.2X. 

1437  *  9HDISS0LVED. 3X . 9HDISS0LVED  ) 

1438  622  FORMAT  (  50X.6H0XYGEN.7X.6H0XYGEN  ) 

1439  625  FORMAT  (  53X.  14HC0NCENTRATI0NS  ) 

1440  630  FORMAT  (  A7 , 2X . F6 . 2 . 7X .F5 . 2 . 7X.F5 . 0. 4X , F5 . 1 . 4X . 

1441  *  F6.2.5X.F6.2  ) 

1442 

1443 


RETURN 

END 


1445 

1446* 

1447* 

1448* 

1449* 

1450* 

1451 

1452 

1453 

1454 
1455* 

1456 

1457 

1458 

1459 

1460 

1461 


SUBROUTINE  ERROR  (  CHECK*  CHAR  ) 

INPUT  DATA  -FILE  ERRORS  FOR  TURBINE  AERATION 
STRUCTURES 


COMMON  /  AA  /  IFILE*  OFILE »  END*  NPOINTS*  NTEST *  DBASE 
COMMON  /  DD  /  QFIRST *  QDEBUG*  QERROR*  QSTOP*  QECHO 
INTEGER  OFILE.  CHECK*  CHAR 
LOGICAL  QERROR 

WRITE  < OFILE*  100)  CHECK.  CHAR 
QERROR  =  .TRUE. 

RETURN 

100  FORMAT  (  '  '  .  '***  INPUT  DATA  ERROR  '  WAS  '*  A4* 

*  2X *  ' **EXPECTING  '»  A4  **'  ) 

END 
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APPENDIX  D:  SAMPLE  INPUT  TO  AND  OUTPUT  FROM 
VENTING  CODE 


CLARKS  HILL  TURBINE  AERATION 

FILE  05  06 

NO  DEBUG 

ENGLISH 

PRINT  OUTPUT 

CENTERLINE  ELEVATION 

AREA  OF  DRAFT  TUBE  EXIT 

INTERVAL  0.1 

BASE  FLOW  TINE 

ALPHA  COEFFICIENT 

BETA  COEFFICIENT 

PRESSURE-TINE  HISTORY 


SEPT  81 


168.0 

672.0 


15.0 

0.33 

0.15 


TINE 


PRESSURE 

15.0 
25.5 
29.1 
25.5 


0.0 
1.4 
2.8 
9.0 

BASE  FLOWRATE 

DATASET  CLARKS  HILL  TURBINE 
TAILUATER  ELEVATION 
TENPERATURE 
AIR  FLOW  49.7 

WATER  FLOW  5280. 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 

TAILWATER  ELEVATION 

TEMPERATURE 

AIR  FLOW  47.0 

WATER  FLOW  5160.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 

TAILWATER  ELEVATION 

TENPERATURE 

AIR  FLOW  39.7 

WATER  FLOW  5193.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 

TAILWATER  ELEVATION 

TEMPERATURE 

AIR  FLOW  29.0 

WATER  FLOW  5190.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 

TAILWATER  ELEVATION 

TEMPERATURE 

AIR  FLOW  20.5 

WATER  FLOW  5193.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 

TAILWATER  ELEVATION 

TEMPERATURE 

AIR  FLOW  10.3 


3700. 
TEST  1 
186.4 
18.9 


TEST  2 
186.4 
18.9 


TEST  3 
186.4 
18.9 


TEST  4 
186.4 
18.9 


TEST  5 
186.4 
18.9 


TEST  6 
186.4 
18.9 


WATER  FLOW  5193.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 
TAILWATER  ELEVATION 
TEMPERATURE 
AIR  FLOW  0.0 

WATER  FLOW  5193.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 


TEST  7 
186.6 
18.9 


TEST  9 


TAILWATER  ELEVATION 

TEMPERATURE 

AIR  FLOW  41.0 


186.8 

18.4 


WATER  FLOW  4887.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE  TEST  10 


TAILWATER  ELEVATION 

TEMPERATURE 

AIR  FLOW  41.0 


186.8 

18.3 


WATER  FLOW  4530. 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE  TEST  13 
TAILWATER  ELEVATION  186.9 

TEMPERATURE  18.3 

AIR  FLOW  39.7 


WATER  FLOW  4144.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE  TEST  14 


TAILWATER  ELEVATION 

TEMPERATURE 

AIR  FLOW  58.1 

WATER  FLOW  3736.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 

TAILWATER  ELEVATION 

TEMPERATURE 

AIR  FLOW  59.7 

WATER  FLOW  3187.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 

TAILWATER  ELEVATION 

TEMPERATURE 

AIR  FLOW  61.1 


185.1 

18.9 


TEST  17 
185.1 
18.9 


TEST  18 
185.1 
18.9 


WATER  FLOW  3511.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE  TEST  21 
TAILWATER  ELEVATION  185.2 

TEMPERATURE  1<M 

AIR  FLOW  53.2 


WATER  FLOW  3951.0 
DO  INITIAL  2.° 

DATASET  CLARKS  HILL  TURBINE  ~EST  22 


TAILWATER  ELEVATION  185.2 

TEMPERATURE  19.0 


r 


TT  rw- 


D4 


AIR  FLOW  49.2 

WATER  FLOW  4378.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 

TAILWATER  ELEVATION 

TEMPERATURE 

AIR  FLOW  50.6 

WATER  FLOW  4747.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 

TAILWATER  ELEVATION 

TEMPERATURE 

AIR  FLOW  53.4 

WATER  FLOW  5046.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 

TAILWATER  ELEVATION 

TEMPERATURE 

AIR  FLOW  48.3 

WATER  FLOW  3420.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 

TAILWATER  ELEVATION 

TEMPERATURE 

AIR  FLOW  39.0 

WATER  FLOW  3431.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 

TAILWATER  ELEVATION 

TEMPERATURE 

AIR  FLOW  24.2 

WATER  FLOW  3429.0 

DO  INITIAL  2.9 

DATASET  CLARKS  HILL  TURBINE 

TAILWATER  ELEVATION 

TEMPERATURE 


TEST  25 
185.2 
19.0 


TEST  26 
185.2 
19.0 


TEST  30 
195.0 
19.0 


TEST  31 
185.0 
19.0 


TEST  32 
185.2 
19.0 


TEST  33 
185.4 
19.0 


AIR  FLOW 
WATER  FLOW 
DO  INITIAL 
STOP 


18.3 

3437.0 

2.9 


